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Abstract 

Conceptual design studies are underway for muon colliders and other high-current 
muon storage rings that have the potential to become the first true "neutrino fac- 
tories" . Muon decays in long straight sections of the storage rings would produce 
precisely characterized beams of electron and muon type neutrinos of unprecedented 
intensity. This article reviews prospects the for these facilities to greatly extend our 
capabilities for neutrino experiments, largely emphasizing the physics of neutrino 
interactions. 

Key words: 

muon colliders, muon storage rings, neutrino factories; PACS13.15.-|-g 



Preprint submitted to Elsevier Preprint 



1 February 2008 



Contents 



1 Overview 



1.1 Introduction 



1.2 Experimental Overview 



1.3 Physics Overview 



2 Deep Inelastic Scattering and QCD Studies 



2.1 Background on Measuring Parton Distribution Functions and QCI^ 



with Non-Polarized Targets 



2.2 Measurement of Quark Parton Distribution Functions 



2.3 Tests of Perturbative QCD| 



2.4 Heavy Quark Production 



2.5 Parton Distribution Functions at Large Bjorken x 



2.6 Examining the Spin Structure of the Nucleon 



2.7 Experimental Setup and Measurement Techniquel 



2.8 Applications of Polarized Parton Distribution Data from z^MCs 

2.9 Studying Nuclear Effects with Neutrinos| 



3 Studies of the CKM Quark Mixing Matrix 



3.1 Introduction 



3.2 Analyses Involving Charm Production: the Extraction of Vcd and 



3.3 Analyses Involving Bottom Production: the Extraction of Vub and 



7 
7 

16 
22 

24 
26 
28 
30 
32 
33 
35 
36 
37 
40 
41 

43 

46 
48 
48 
49 
49 



3.4 Vcd/^cs via Diffractive Charmed Vector Meson Production 



3.5 Improved Knowledge of the CKM Matrix from uMCs 



Precision Electroweak Studies 



4.1 Introduction 



4.2 Elastic and Quasi-Elastic Neutrino-Electron scattering 



51 



2 



43 sin'^ 9w from Deep Inelastic Scattering 



4~4 Summary on sin"^ 9w Measurements at vM.Cs 



5 Rare and Exotic Processes 



5.1 New Physics Sensitivity 



5.2 Studies of Low Energy QCD 



5.3 Conclusions on Rare and Exotic Processes at i/MCs 



6 Charm Decay Physics 



6.1 Introduction 



6.2 Theoretical Motivation for Charm Physicsl 



6.3 Probing Strong Interactions through Charm Decays 



6.4 Searches for New Physics in Charm Decays 



6.5 Summary on Charm Decay Physics at i^MCs 



7 Neutrino Oscillation Experiments with a Muon Storage 
Ring/Neutrino Factory 



7.1 Status of Neutrino Oscillations at the Time of z^MCs 



7.2 Oscillation Experiments at t/MCs| 



7.3 Matter Effects 



7.4 Detector Considerations 



7.5 Conclusions on Neutrino Oscillation Studies at vMCs 



8 Summary 



9 Acknowledgments 



References 



60 
68 
69 
69 
74 
78 
78 
78 
79 
80 
84 
90 

90 
90 
91 
94 
95 



99 
99 
100 



3 



List of Figures 



1 Example neutrino event spectra for and Ue from a 50 GeV 



negative muon beam from a neutrino factory 10 


Solid curves 


indicate the spectra for decays at zero degrees in the center 


ot mass system, 'i'his is the spectrum expected for a detector 


located very far from the muon decay region. The dashec 


1 


curves indicate the spectra for decays within the forward 




hemisphere in the center of mass frame. This is what would be 


expected for a detector close enough to the muon decay region 


to subtend an angle of I/7. 






2 Example of a general purpose neutrino detector 


6 . Its scale 


is illustrated by a human figure in the lower left corner. The 


neutrino target is the small horizontal cylinder at mid-height 


on the right hand side of the detector. Its radial extent 




corresponds roughly to the radial spread of the neutrino 


pencil beam, which is incident from the right hand side. The 


illustration is partially schematic in that the geometries c 


)1 


the calorimeters and dipole magnet have been simplified for 



illustrative purposes. 



3 Conceptual illustration of the vertex tagging superiority 



expected at i^MCs over that with collider experiment 




geometries. (The figure is reproduced from 


14 


, which used 




the terminology "IVIUKIINE", for iVlUon King 


Neutrino 




Experiment, instead of z/MC.) Neutrino targets could have 




a vertex plane ol CCD pixel detectors spaced at intervals 




approximately one millimeter. For comparison, 


the VXD3 




vertexing detector at the SLD experiment at ISLACJ, generally 


regarded as the best existing vertex detector in a collider 


experiment, has its two innermost CCD tracking planes at 2.8 


cm and 3.8 cm from the interaction point (IP) 


A schematic 


of a one-prong decay has been drawn to illustrate the 




advantages ot cfoscly spaced vertex detectors, i'br clarity ol 




illustration, the kink defiection angle has been 


drawn much 




larger than would be typical, 'i'he 2 cm dist 


ance to decay for 


the charmed meson corresponds to the 


average boosted 





lifetime for a 120 GeV D+. 



4 Fractions of the total neutrino-nucleon cross section involving 



production of heavy flavors in the final state, for an isoscalai 



target and as a function of neutrino energy|8|. 'I'he plotted 



production fractions are for charged current charm production 



[RcE, and RdoE is the leading order approximation), neutral 



current production of a charm- ant icharm pair (Rcce), charged 



current B production from a u quark [RbuE, again with 



RbuioE as the leading order approximation), charged current 



B production from a c quark {RbcE) and neutral current 



production of a bh pair [RbbE) 



Heavy quark fractional rates vs neutrino energies corresponding 
to l:''ig. 4, but for antineutrinos rather than neutrinos 8 . 





6 Model-dependent decomposition of nucleon 


spin into 




contributions from quarks and gluons. 'i'he two vertical dotted 




lines show the naive QPM expectation and the results 


from a 


NLO fit to most of the available data on g\ . 





7 Simulated predictions for the ratios of both the F2 and xh^ 



structure functions between calcium and deuterium, plotted as 



a function of Bjorken x and assuming the theoretical model ol 
Kulagin |45| and data samples from a z/MC. The sizes of the 



error bars correspond approximately to 1-year exposures foi 



both targets at the 50 GeV t/MC of Table 1 



8 Feynman diagrams contributing to the v — e scattering 



processes of Eqs. 62 to 66: (a) NC v-e elastic scattering (for 



Eqs. 62, 64, 65 and 66) (b) CC u^-e scattering (for Eq. 65), (c 



UC annihilation (for equation 66), and (d) inverse muon 



decay (for Eq. 64)^ 





9 Integral cross section for neutrino-electron scattering processes 


above any chosen cut on the inelasticity variable, y > r/cut, and 


assuming = 30 GeV. 







10 Ghange in the integral cross section with respect to sin 6w for 


neutrino-electron scattering processes above y > ?/cut, assuming 


K = 30 GeV. 







11 Statistical uncertainty in sin Ow for any chosen value of ?/cut, 


from neutrino-electron scattering in beams of either Vf., 


u 


a 




or z/g. An integrated luminosity of 10"*° cm ^ at a beam energ 


y 


of Ey = 30 GeV has been assumed. The values can be scalec 




to other neutrino energies by noting that the measurement's 




statistical uncertainty, for a given integrated neutrino flux 








through a specified detector, is proportional to the inverse 






square root of the average energy.] 









12 Same as Fig. 11 except the event statistics has been summed 


over the experimentally indistinguishable contributions from 


the two neutrino components in the //"""-induced i/^z/g beams 




and the fi -induced i/^z/g beams produced from muon storage 


rmgs. 





13 Monte-carlo generated simulation of a high energy neutrino- 
electron scattering interaction in liquid helium. 'I'he distorted 
scale is indicated by the 10:1 ratio in the relative lengths ol 
the X and z axes. The solid lines are electrons and positrons| 
and the dashed lines are photons that would not be seen in the 

detector. The view is perpendicular to both the beam direction 

and to a 0.1 'lesla magnetic held. See text for further details" 





14 Expected future restrictions in the size of the allowed regions 


in a plot of W mass vs. top quark mass, following experiments 


at 'ieV 33 and at a z/iVl(J, as illustrated by the thicknesses 




of the bands and ellipse. The exact positions of the shadec 


regions are for illustrative purposes only. 







15 Plot 


s for limits 


on th 


e (from top to bottom) coupling values 






UeL 


'2 


U^LUeL 


and 


UaL ^ as a function of Lq mass for one 






year 


ol 


runnmg 


under the assumptions given in the text, 'i'he 




solid, d 


ashed, and dotted lines represent uMC energies of 10, 




20, and 250 GcV respectively, 'llic plots on the lelt show limits 




for L = 


12 m, and the plots on the right show limits for L = 


1 


km. 





16 'i'he upper plot shows and z/^ event rates per GeV as c 




function of E/L for two different experiments. The lower 




plot shows the L/E dependence of the oscillation probability. 


assuming the largest dm"^ is 2,3.5,or 5 xlO~'^ eV^/c"'. 



1 Overview 



6 



1.1 Introduction 



Muon colliders have been proposed to provide lepton-lepton collisions while 
circumventing the energy limitations on electron-positron storage rings caused 

by synchrotron radiation. The larger muon mass suppresses synchrotron ra- 
diation energy losses by a factor m^/nv"^ ~ 5 x 10^"*^° relative to those of a 
circulating electron beam of the same energy and, incidentally, also opens up 
promising possibilities for s— channel Higgs boson production[l]. 

Recent feasibility and design studies for future muon colliders [2.1] have begun 
to focus attention on the exciting physics possibilities for experiments using 
neutrino beams from the decays of the circulating high energy muons. This 
report explores the potential for a "neutrino experiment at a muon collider" , 
or z/MC for short. A i/MC program could operate either parasitically during 
a colliding beam experiment; or it could be installed as part of a program in 
neutrino physics at a dedicated muon storage ring. 

Amongst the potential physics topics for i/MCs, neutrino oscillations have 
garnered the most intense experimental and theoretical activity, and particular 
possibilities for long baseline oscillation experiments exploiting a muon storage 
ring are covered elsewhere [3, 4]. In this report, we wish to also highlight the 
superb capabilities of neutrinos as probes of the strong and weak interaction 
dynamics of quarks and the parton structure of nucleons, as well as the power 
of a vMC in searches for evidence of new types of weak interactions. 

The remainder of this section lays out the expected experimental parameters 
and capabilities of a vMC and provides concise overviews for the more detailed 
physics discussions that follow. 

1.2 Experimental Overview 

1.2.1 High Current Muon Storage Rings 

Recent ideas for neutrino experiments at either muon colliders [5,6] or ded- 
icated neutrino factories [7] represent reincarnations of earlier proposals for 
neutrino experiments at muon storage rings that date back at least to the 
1960 's. The essential advantages of modern i/MCs derive from the very large 
muon currents that might be produced and stored using the technologies devel- 
oped for muon colliders. Current design scenarios for muon colliders [2,1] and 
neutrino factories [10,4] envision of order 10^^ positive and negative muons 
per year circulating and decaying in the storage ring. 

Neutrinos from decays in the longest straight sections of the storage ring will 
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description 



z^— factory Higgs-factory top-factory 



muon energy, 



.20 



/i±/ycar 

flight time to beam dump ti)/7r^ 
ring circumference, C[m] 
straight section (SS) length Zss[m] 
fractional SS length, [fss= hs/C] 
/x+/year in SS, [iV-= /..iV^/lO^o" 
u from SS/ycar /lO^*^ 
angular divergence [/7 • (56* j^] 
1/ angular divergence [mrad] 
N^^ [events /yr /g/ cm^ ] 
target thickness[g cm~^] for 10^°events 
N^^ [events /yr /kT / ( lO^km) ^] 



20 GeV 

3.0 
no dump 
300 
90 
0.30 
0.90 
1.8 
1 

5.3 
3.8 X 10^ 

2600 
1.2 X 10^ 



50 GeV 

6.0 
no dump 
345 
40 
0.12 
0.72 
1.4 
1 

2.1 
6.5 X 10^ 

1500 
1.4 X 10^ 



175 GeV 

6.0 
no dump 

900 

110 

0.12 

0.72 

1.4 
1 

0.60 
2.7 X 10^ 

370 
6.2 X 10^ 



description 



muon energy, 

/x±/year [a^^/IO^^ 
flight time to beam dump to/jT^ 
ring circumference, C[m] 
straight section (SS) length Zss[m] 
fractional SS length [fss= hs/C] 
//+/year in SS [n^^= f.sN^/lO^'^ 
V from SS/year [/lO^" 
V angular divergence [/7 • SQ^ 
V angular divergence [mrad] 
N^^ [events Iy^IzI ] 
target thickness[g cm~^] for 10^°events 
N^^ [events /yr /kT/ ( lO^km) ^] 



frontier 2nd generation 



500 GeV 


5 TeV 


3.2 


3.6 


0.5 


no dump 


2000 


15 000 


150 


450 


0.12 


0.03 


0.38 


0.11 


0.30 


0.22 


10 


1 


2.1 


0.021 


2.3 X 10'^ 


1.0 X 10^ 


430 


100 


5.0 X 10^ 


2.2 X 10^° 



Table 1 

Neutrino fluxes and event rates for representative example parameter sets for ded- 
icated neutrino factories or muon colliders [9] spanning the energy range = 20 
GeV to 5 TeV. The angular divergence scaling factor, 7 • bQ^ is the factor by which 
the divergence of the parent muon beam increases the neutrino beam's angular di- 
vergence beyond the characteristic size, = I/7, expected for a divergenceless 
muon beam. 



Neutrino Fluxes from a 50 GeVMuon Beam 




neutrino energy, GeV 

Fig. 1. Example neutrino event spectra for and z^e from a 50 GeV negative muon 
beam from a neutrino factory[10]. Solid curves indicate the spectra for decays at 
zero degrees in the center of mass system. This is the spectrum expected for a 
detector located very far from the muon decay region. The dashed curves indicate 
the spectra for decays within the forward hemisphere in the center of mass frame. 
This is what would be expected for a detector close enough to the muon decay 
region to subtend an angle of I/7. 

emerge in intense coUinear beams that are highly suitable for experiments. 
Beams from such production straight sections should provide many orders of 
magnitude higher event rates than considered in the early versions of muon 
storage rings and, indeed, should be considerably more intense than today's 
"conventional" neutrino beams produced from tt/K decays. No need exists 
for a conventional beam's muon shielding berm, and detectors can be placed 
relatively close to the end of the production straight section. Coupled with the 
relativistic kinematics of muon decay, this permits the possibility of detectors 
only tens of centimeters across and allows for the use of high precision vertex 
detectors as active neutrino targets. 
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target purpose 


general 


polarized 


V — e 


material 




solid tl U 


hquiQ 0/24 


mean density [g/cm*^] 


U.O 


U.ZD i 


U. / i / 


length [m] 


o 
Z 


0.5 


on 
zO 


thickness [g cm 


iUU 


1 Q /I 

io.4 


1 /I Qn 


radius [cm] 


on 


on 


on 
zU 


mass [kg] 


iZD 


iu.o 


1 Qnn 
ioUU 


integrated luminosity [fb ^] 


c n V 1 n6 

D.U X iU 


o.i X iU 


o.D X iU 


DIS events/year at 50 GeV 


1 1 ^/ A n8 
i .{ X iU 


i.U X iU 


11.^1 nlO 
i.i X iU 


DIS events/year at 175 GeV 


2.7 X 10^ 


3.6 X 10^ 


3.8 X 10^° 


z/e events/year at 50 GeV 


2 X 10^ 


NA 


3 X 10*^ 


Ve events/year at 175 GeV 


7 X 10^ 


NA 


1 X lO'^ 



Table 2 

Specifications, integrated luminosities and event rates for the high rate neutrino 
targets discussed in this report, assuming the 50 GeV and 175 GeV muon storage 
ring parameters of Table 1. The target is assumed to be situated 100 m (350 m) 
downstream from the center of the 50 GeV (175 GeV) production straight section. 

Additional physics advantages over tt/K decay neutrino beams will result from 
the unique and precisely understood flux composition of the uMC beams. 
Negative and positive muons decay according to 

fi~ ^ i^^, + i^e + e", 

fx'^ ^ Uf, + i^e + e+, (1) 

producing pure 2-component neutrino beamsQ via the perhaps best under- 
stood of all weak decay processes. These beams will be designated as z/^z/e or 
Ufj^Ue, respectively, in the rest of this report. 

Experimental requirements for the two broad classes of neutrino physics at 
z/MCs differ greatly, chiefly because the experiments would be conducted at 
very different baseline distances from beam production to the detector. Ex- 
periments for neutrino interaction physics will be conducted as close to the 
muon ring as possible ( "short baseline" ) in order to maximize event rates and 
to subtend the neutrino beam with a target of small transverse cross sec- 
tion. On the other hand, the choice of baseline for neutrino oscillation studies 

^ We implicitly assume here the absence of a significant lepton family number 
violating decay of the type /i~ e^z^g^^^t but caution that the current experimental 
limit on the branching fraction for this decay is only 1.5%. This limit will clearly 
be greatly improved upon from the consistency of the observed z^MC spectra with 
predictions. 



10 



will be dictated by the specific range of possible oscillation parameters un- 
der investigation, as discussed further in Chapter 7. Oscillation parameters of 
current interest motivate the use of very long baselines, even extending to the 
possibility of transcontinental experiments [7]. 

As an important caveat on the contents of this report, it should always be 
borne in mind that the ambitious technologies of these high current muon 
storage rings still only exist at the feasibility or early design study stage and 
it is by no means guaranteed that realizable devices will appear anytime soon. 
Nevertheless, recent progress has been impressive, and the pace of R&D is 
accelerating. The reader is referred to the specialist literature for a more thor- 
ough overview of the technological challenges in building a muon collider[2,l] 
or neutrino factory[10,4]. 



1.2.2 Event Rates 



Event rates in all uMC experiments will be dominated by the charged current 
(CC) and neutral current (NC) deep inelastic scattering (DIS) of neutrinos or 
antineutrinos with nucleons {N — p or n): 



ut + N ^ r + X 

vt + N + X 



(i/,(i/,)-NC), 



(2) 



where i = e or fj, and X represents a typically multi-particle hadronic final 
state. Neutrino-nucleon DIS cross sections scale with neutrino energy E^, to 
a good approximation for neutrino energies above a few GeV, with numerical 
values of [12]: 



cc] 






NC 




2.1 


CC 




3.4 


NC J 




^1.3 J 



X E^[GqV] fb. 



(3) 



At the many-GeV energies of i^MCs, vN DIS is well described as the quasi- 
elastic scattering of neutrinos off one of the many quarks or antiquarks inside 
the nucleon through the exchange of a virtual W or Z boson: 

M^e) M^e) + Q (NC), (4) 

i^e + Qd (qu) ^r + q'^ {q',) {u - CC) , (5) 

i>i + qu{qd)^i^ + q'AQ'u) (V-CC). (6) 
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All quarks q participate in the NC process. The CC interactions change quark 
flavor, with neutrino interactions producing w-type and d-type flnal state 
quarks, q!^andq'^, from d-type and M-type targets, q^ and Antineutrinos 
participate in the charge-conjugate processes. Much of the richness of neu- 
trino interaction physics derives from the variety of processes contained in 
Eqs. 4-6. 



1.2.3 Neutrino Production Spectra and Event Rates in Detectors 

Neutrino flux spectra at uMCs will be precisely predictable since the decay of 
muons is a well-understood purely electroweak process. Characteristics of the 
parent muon beam in the production straight section can be reliably calculated 
and modeled through a knowledge of the focusing magnet lattice and through 
beam monitoring. Calibration of the muon energies in the storage ring might 
reach the level of a few parts per million fractional uncertainty [11]. 

Due to the differing angular coverages, the neutrino spectrum seen by an oscil- 
lation detector at a long baseline will differ from that seen by detectors placed 
at short baselines to study interaction physics. Long baseline detectors will 
sample the very forward-going neutrinos, at angles in the muon rest frame(6'') 
and laboratory frame(6') close to 6*' = 6' = 0, while detectors close to the pro- 
duction straight section will instead accept a production solid angle bite that 
is comparable to the boosted forward hemisphere of the decaying neutrinos, 

U 2^y,_smy, i/7 -£;^[TeV]- ^> 



Figure 1 gives an illustrative example of the neutrino spectra at z/MCs for 
detectors at both short and long baselines[10], and Table 1 gives beam and 
event rate parameters for several other vMC scenarios. Further explanation 
for the choices of storage ring parameters in Table 1 and a derivation for the 
following simple numerical expressions for event rates used to fill Table 2 are 
provided elsewhere [9]. 

For short baseline detectors. 



N'^ [events/yr/g/cm-2] = 2.1 x 10"^^ x E^[GeV] x Nl'[yr-\ (8) 

where A^*'' is the number of neutrino interactions per year of running per 
g-cm~^ of a cylindrically symmetric target centered on the beam, is the 
muon beam energy and is the number of forward-going muons (as op- 
posed to muons circulating in the opposite direction in, e.g., a collider ring) 
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decaying in the production straight section per year. Equation 8 assumes the 
parent muon beam to have an angular divergence in the production straight 
section that is small compared to the 50 y = I/7 natural divergence of the neu- 
trino beam. This will normally be the case [13] unless the choice of straight 
section is in the final focus region of a coUider storage ring. Table 2 accounts 
simplistically for this exception (in this case for the illustrative parameter set 
at 500 GeV) by increasing the angular divergence of the neutrino beam by a 
simple scale factor. In this circumstance, the angular coverage of the target 
would need to be increased by this same scale factor in order to retain the 
event rate predicted by the parameter A^*^. 

In contrast to short baseline detectors, the event rate in long baseline detectors 
is not sensitive to the geometry of the detector since the entire detector will 
always be bathed uniformly by the forward-going neutrino flux. The number 
of interactions in the detector will vary in proportion to the target mass M 
and inversely as the square of the baseline length L. This leads to a definition, 
analogous to Eq. 8, for the event rate benchmark N'''^: 



-20 



events/yr/kT/ ^10^ km^ 



1.6 X 10-2" xA^r[y,--i]x(£,[G'ry])= 



(7 • sey 



(9) 



where N^^ is the number of neutrino interactions per kiloton-year of running 
with a target centered on the beam at a 1000 km distance from the production 
point. The previously discussed angular divergence scaling factor, 7 • 59, has 
been explicitly included. 

The event rates given in Table 2 are truly impressive. Samples of thousands of 
events per kiloton might be recorded at oscillation experiments with baselines 
as long as thousands of kilometers. For neutrino interaction physics, samples 
as large as 10 billion events can be reasonably contemplated in compact targets 
close to the production straight section. Equation 7 shows that the radial 
extent of such targets can be as small 10-20 cm. 



1.2.4 Detector Design Considerations for vMCs 

Event rates for oscillation experiments will probably be less of an extrapola- 
tion from today's experiments than will be the case for interaction experiments 
due to the compensating rate decrease at the expected longer baselines. Cor- 
respondingly, the innovations in neutrino detector design required to upgrade 
to the neutrino beams at uMCs are likely to be rather less substantial for oscil- 
lation experiments at long baselines than for interaction physics experiments. 

Two significant changes expected for the design of oscillation detectors for 
i/MCs are that (i) the 2-component beams provide strong motivation for a 
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Fig. 2. Example of a general purpose neutrino detector [6]. Its scale is illustrated by 
a human figure in the lower left corner. The neutrino target is the small horizontal 
cylinder at mid-height on the right hand side of the detector. Its radial extent 
corresponds roughly to the radial spread of the neutrino pencil beam, which is 
incident from the right hand side. The illustration is partially schematic in that the 
geometries of the calorimeters and dipole magnet have been simplified for illustrative 
purposes. 

magnetic spectrometer to distinguish muon charge signs; and (ii) larger de- 
tector masses might be financially justified in order to fully exploit the large 
financial investment in the muon storage ring. Design considerations for de- 
tectors for oscillation z/MCs are discussed in more detail in Chapter 7. 

In contrast to oscillation experiments, the increase in neutrino yield for uMCs 
relative to beams from pion decays as well as the coUimation of the neutrino 
beams will allow the use of compact, speciahzed targets surrounded by high 
performance detectors. These detectors must operate at high rate in order to 
cope with the data sets implied in Table 1. Considerable thought must be 
given to triggering, data acquisition, event reconstruction and data handling 
considerations. 

Figure 2 provides an example[6,9] of the sort of high rate general purpose neu- 
trino detector that would be well matched to the intense neutrino beams at 
uMCs. The neutrino target is one meter long stack of CCD tracking planes 
represented by the small horizontal cylinder at mid-height on the right-hand 
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interaction 




,00' 



00 000 0000 0000000 000° 

Oq^^ 1-prong D+ decay (cartoon) 

A) vertexing at MURINE 



B) best collider vertexing so far: (SLD detector, same scale) 



IP 



Fig. 3. Conceptual illustration of the vertex tagging superiority expected at z^MCs 
over that with collider experiment geometries. (The figure is reproduced from [14], 
which used the terminology "MURINE", for MUon Ring Neutrino Experiment, 
instead of vMC) Neutrino targets could have a vertex plane of CCD pixel detectors 
spaced at intervals of approximately one millimeter. For comparison, the VXDS 
vertexing detector at the SLD experiment at SLAC, generally regarded as the best 
existing vertex detector in a collider experiment, has its two innermost CCD tracking 
planes at 2.8 cm and 3.8 cm from the interaction point (IP). A schematic of a 
one-prong decay has been drawn to illustrate the advantages of closely spaced 
vertex detectors. For clarity of illustration, the kink deflection angle has been drawn 
much larger than would be typical. The 2 cm distance to decay for the D'^ charmed 
meson corresponds to the average boosted lifetime for a 120 GeV D~^. 

side of the detector in Fig. 2. Its 10 cm radial extent could correspond to, 
e.g., the 0.2 mrad divergence of the neutrino beam originating from a 500 GeV 
muon beam 500 m upstream of the target. The scale of the entire detector is 
illustrated by a human figure in the lower left corner, emphasizing the strik- 
ing contrast in target size with the kiloton-scale coarse-sampling calorimetric 
targets often used for past and present high rate neutrino experiments. 

The CCD target in Fig. 2 contains 750 planes of 300 micron thick silicon 
CCD's, corresponding to a mass per unit area of approximately 50 g-cm~^; 
this translates to 2.5 radiation lengths or 0.5 interaction lengths. Scaling to 
different target lengths and radii should be straightforward without altering 
the basic design of the surrounding detector. 
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Besides providing the mass for neutrino interactions, the tracking target al- 
lows for precise reconstruction of the event topologies from charged tracks, 
including event-by-event vertex tagging and reconstruction of those interac- 
tions containing heavy flavor final states. The fixed target geometry of i/MC 
vertex detectors allows for much more frequent sampling than is possible in 
collider detectors: Fig. 3 gives a schematic comparison between the charm 
vertexing capabilities of the CCD detector of Fig. 2 and the current best ver- 
texing detector in a collider experiment [14]. 

The CCD target is backed by a hermetic detector reminiscent of many col- 
lider detector designs. An enveloping time projection chamber (TPC) pro- 
vides track- following, momentum measurements, and particle identification 
for essentially all charged tracks emanating from the interactions. Optionally, 
further particle identification might be available from a mirror that refiects 
Cherenkov light to an instrumented back-plane directly upstream from the 
target. The mirror is backed by electromagnetic and hadronic calorimeters 
and, lastly, by iron-core toroidal magnets for muon identification. 

Other possible specialized high rate neutrino target and detector possibilities 
include polarized solid protium-deuterium targets for spin physics (Sec. 2.6) 
and nuclear targets (Sec. 2.9) for studies of A dependence. A more massive 
tracking liquid target (Sec. 4) would be suitable for precision electroweak 
physics using neutrino-electron scattering. Table 2 provides a summary of 
some of the characteristics for examples of each of the three high rate target 
types discussed in this section and also gives plausible but very approximate 
integrated luminosities and event sample sizes for of the illustrative 50 GeV 
and 175 GeV beam parameters in Table 1. 

1.3 Physics Overview 

This overview motivates and introduces the more detailed discussions that 
follow on: deep inelastic scattering and quantum chromo dynamics (Sec. 2), 
quark mixing (Sec. 3), precision electroweak tests (Sec. 4), rare and exotic 
processes (Sec. 5), charm physics (Sec. 6) and neutrino oscillations (Sec. 7). 

Before proceeding, we note that much of the interesting physics involves as- 
pects of CC and NC charm and beauty production (and hence the motivation 
for active vertex detectors as targets). 

As well as a CKM physics program that complements those from B and K 

factories and from precision W boson branching fraction measurements at 
colliders, b and c production at a z/MC allows precisions tests of QCD near 
heavy fiavor thresholds, permits sensitive probes for new physics such as fiavor- 
changing neutral currents, and provides a novel, very high statistics sample of 



16 



charmed hadrons. 

Figures 4 and 5 show heavy quark production fractions and indicate that, 
given the expected multi-biUion inclusive event samples, very high statistics 
can indeed be accumulated for both c and b final states at sufficiently above 
the relevant energy thresholds. 

1.3.1 QCD and Deep Inelastic Scattering 

Historically, neutrino experiments have made major contributions to our un- 
derstanding and verification of the both the QCD theory of strong interactions 
and the constituent components of protons and neutrons. The extrapolation 
of present experimental statistics consisting of 10^~^ events to the expected 
10^~^° well reconstructed DIS events at z/MCs might well provide the best ever 
experimental laboratory for studying QCD and the structure of the nucleon 
through a scattering process. 

Both traditional and novel areas for potential study will be discussed in Sec. 2. 
They include: (1) one of the most precise and theoretically sound measure- 
ments of the strong coupling constant, a,^; (2) stringent consistency checks 
for the predictions of perturbative QCD; (3) detailed flavor and spin depen- 
dence nucleon structure functions using both CC and NC probes; (4) precise 
tests of QCD near the c and h quark heavy flavor transitions; and (5) the first 
systematic studies of QCD in nuclear environments probed by neutrino and 
antineutrino beams. 

1.3.2 The CKM Quark Mixing Matrix 

Some of the most important high-rate measurements will involve the Cabibbo- 
Kobayashi-Maskawa (CKM) mixing matrix that characterizes CC weak inter- 
actions of quarks. This topic is discussed in detail in Sec. 3. 

Neutrino-nucleon DIS offers unique and systematically independent measure- 
ments of CKM matrix elements since it uses a high virtual W probe cou- 
pling directly to quarks rather than relying on the complex interplay of weak 
and strong interactions that is inherent in hadron decay. For sufficiently high 
energy and event rates, four of the nine CKM matrix element amplitudes - 
|Ks|, iKtbl and \Vcb\ - are directly probed through c and h production. 

The higher momentum transfers from the external W probe allow for a cleaner 

theoretical interpretation that requires only relatively small corrections from 
perturbative QCD. As a further theoretical advantage, the measurements are 
semi-inclusive - i.e. summing over all final states with single charm or beauty 
production - and thus do not suffer from uncertainties in hadronic branching 
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Fig. 4. Fractions of the total neutrino-nucleon cross section involving production 
of heavy flavors in the final state, for an isoscalar target and as a function of neu- 
trino energy [8]. The plotted production fractions are for charged current charm 
production {RcE, and RdoE is the leading order approximation), neutral current 
production of a charm-anticharm pair (Rcce), charged current B production from 
a u quark {RbuE, again with RbuioE as the leading order approximation), charged 
current B production from a c quark {RtcE) and neutral current production of a bb 
pair (RbbE)- 

ratios. 

The fractional production rates shown in figures 4 and 5 for i>—CC c and b 
production at high energies are of order \Vcdf and respectively, where 

(Vcd, Vub) are the {d ^ c, u b) Cabibbo-Kobayashi-Maskawa (CKM) matrix 
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Fig. 5. Heavy quark fractional rates vs neutrino energies corresponding to Fig. 4, 
but for antineutrinos rather than neutrinos [81. 



elements. Other significant contributions to heavy flavor production are pro- 
portional to \Vcsf and \Vcbf, where {VcsyVcb) are the (s — >• c, c — 6) CKM 
elements. 



The relatively clean theoretical interpretations and large samples of flavor- 
tagged events, particularly for charm production, should allow impressive mea- 
surements of the absolute squares for several of the elements in the CKM quark 
mixing matrix. Estimated precisions in determining the CKM matrix elements 
are summarized in Table 3. 
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U.94(5 


r\ r\AO 

U.U4o 


1.45 X lU 




±0.16% 


±2.1% 


±60% ^ 0(3%) 


c 


0.050 


1.08 


1.6 X 10-3 




±14% ^ 0(1%) 


±31% ^ 0(3%) 


±10.5% ^ 0(3%) 



Table 3 

Absolute squares of the elements in the first two rows of the Cabbibo-Kobayashi- 
Maskawa (CKM) quark mixing matrix, along with their uncertainties when no uni- 
tarity constraints are applied [15]. The second row of the entry for each element 
gives current percentage one-sigma uncertainties in the absolute squares and pro- 
jections for the uncertainties after analyses from a vMC operating with neutrino 
energies well above the B production threshold. 

Perhaps the most interesting potential measurement outlined in Sec. 3 is the 
determination of \Vub\ to better than 5%, perhaps eventually reaching 1%. This 
is an order of magnitude better than the current uncertainty and might well be 
better than will be achieved in any other single measurement at, for example, 
a B factory. B production should also allow for an extraction of \Vcb\ at the 
few percent level that is systematically different, in both its experimental and 
theoretical aspects, from studies of decay processes and that is comparably 
accurate to the anticipated future measurements using decays. 

The matrix clement \Vcd\ is already best measured from CC charm production 
in today's neutrino experiments, based on event samples of several thousands 
of events. The present accuracy is mainly limited by statistics and uncertain- 
ties in charmed hadron production and decay characteristics. It is clear that 
the accuracy in \Vcd\ would be vastly improved from the analysis of hundreds 
of millions of vertex-tagged charm events in a high performance detector. 

1.3.3 Precision Electroweak Physics 

Section 4 demonstrates that z/MCs should be able to provide two types of preci- 
sion measurements of the weak mixing angle sin^ 9]^: from the ratio of neutral 
current (NC) to charged current (CC) DIS events and also from neutrino- 
electron scattering. 

Both types of determinations require a large extrapolation in event statistics 
and experimental technique from today's best neutrino results. They will allow 
vigorous consistency checks of the Standard Model and provide sensitivity to 
several potential possibilities for new physics. 
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With its huge statistics, the DIS measurement of sin^ 9w will eventually be sys- 
tematically limited by theoretical hadronic uncertainties but it should anyway 
become several times more precise than today's best neutrino measurements, 
which are already equivalent to about a 100 MeV uncertainty on the W mass. 

By contrast, no significant theoretical uncertainties enter into neutrino-electron 
scattering - a simple scattering process between two elementary point particles 
- and so the measurements will be limited only by statistics and experimental 
ingenuity. One can contemplate neutrino-electron scattering event samples 
as large as 10^ events using a dedicated detector with parameters like those 
given in Table 2. This would correspond to impressive statistical uncertainties 
in sin^^^ of order 10^^ and sensitivity to new contact interactions at energy 
scales up to approximately 25 TeV. The biggest experimental challenges may 
come from normalizing the neutrino beam flux. If the experimental uncer- 
tainties could be reduced to the extremely challenging level of the statistical 
uncertainties then this process holds the potential for measurements of sin^ 6w 
that might potentially be as good as or better than the best current measure- 
ments from collider experiments. 

1.3.4 Rare and Exotic Processes 

A i/MC provides a facility for several unique searches for physics beyond the 
Standard Model and provides a venue for observing a number of very rare 
Standard Model processes, as outhned in Sec. 5. Examples of the former in- 
clude flavor changing u ^ c, d ^ b, and s ^ b neutral currents and isosinglet 
electron and muon-type neutral heavy leptons. A list of the latter includes 
annihilation and the scattering of virtual W and Z bosons from quasi- 
real photons in a nuclear Coulomb held; these serve as weak analogs of e+e" 
and 77 physics. Sensitivity to high energy scales through mixing of or Z 
with higher mass propagators is limited due to the characteristically weak de- 
pendence of the new propagator mass reach as the fourth root of experimental 
statistics. 



1.3.5 Charm Decays 

A uMC should function as an efficient factory for the study of charm decays, 
with a clean, well reconstructed sample of several times 10^ charmed hadrons 
produced in 10^° neutrino interactions. Section 6 points out several interesting 
physics motivations for charm studies at a uMC. Measurement of charm decay 
branching ratios and lifetimes are useful both for QCD studies and for the 
theoretical calibration of the physics analyses on B hadrons. Charm decays 
also provide a clean laboratory to search for exotic physics contributions since 
the Standard Model predicts tiny branching fractions for rare decays, small 
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CP asymmetries and slow D'^ — > D° oscillations. 

The charge of the final state lepton in CC- induced charm production from neu- 
trinos tags charm quarks vs. antiquaries with high efficiency and purity. This 
tag is of particular benefit to oscillation and CP studies, as is the expected 
precise vertexing reconstruction for the proper lifetime of decays. Section 6 
shows it to be quite plausible that a uMC could provide the first observations 
of both D^ — D^ mixing and CP violation in the charm sector and additionally 
provide some context for their proper interpretation. 

1.3.6 Neutrino Oscillations 

The potential for long basehne oscillation experiments is the uMC topic that 
is currently of most interest to the high energy physics community. However, 
both the experimental and theoretical status of neutrino oscillations are in 
such a state of fiux that long-term predictions for i/MCs can be stated only in 
general terms. 

Long baseline experiments at z/MCs might well provide a definitive follow-up to 
the recent intriguing evidence for neutrino oscillations. There will presumably 
already have been some progress in the verification or refutation of today's 
oscillation signals by the time long baseline uMCs come on-line. Even so, 
i/MCs will still clearly be important for more probative follow-up studies to 
characterize the form and phenomenology of any observed oscillation signals. 
For example, z/MCs might help to determine whether an observed oscillation 
signal is consistent with mixing between three neutrino families or whether 
a fourth, sterile neutrino is required, as is discussed in Chapter 7. In the 
former case, long baseline i/MCs can nail down the values of the mass-squared 
differences, search for CP violation and look for matter effects in oscillations in 
the Earth's interior. Clearly, the spectrum of possible studies would be richer 
still in the case of 4-neutrino mixing. 



2 Deep Inelastic Scattering and QCD Studies 

Starting from SLAC electron scattering experiments in the late 60's, proceed- 
ing through the CERN and Fermilab neutrino and muon experiments of the 
70's through 90's, and continuing with HERA experiments still underway, 
deep inelastic scattering (DIS) has provided us with an increasingly accurate 
picture of the partonic structure of the nucleon. Moreover, DIS has served and 
still serves as one of the best test-beds for perturbative QCD. 

A uMC could take the physics of DIS to a new level by: (1) providing the 
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statistical power to extract all six structure functions for v and V beams on 
proton and deuterium targets; (2) allowing for low mass, high acceptance spec- 
trometers with vastly improved resolution over present calorimetric detectors; 
(3) creating naturally redundant measurements through simultaneous mea- 
surement of electron and muon final scattering final states; (4) generating the 
rate and small beam spot size required for the first polarized neutrino tar- 
gets; (5) permitting use of active vertexing targets for systematic studies of 
heavy flavor production, and (6) facilitating the use of a large array of nuclear 
targets. 

Examples of physics topics that would emerge include: 

• Definitive proton parton distribution functions (PDF) for x > 0.01. The 
X > 0.01 behavior of parton densities at accessible uMC controls the 
cross section behavior of the highest energy scale physics of the Tevatron 
and LHC Understanding subtle deviations caused by new physics requires 
precise control of PDF systematics. A z/MC will have the statistical power 
and the systematic redundancy checks to generate a complete PDF set from 
a single experiment. 

• A test of QCD to next-to-next to leading order (NNLO). Few experiments 
quantitatively test QCD beyond leading order in the coupling ag (Q^). It 
is frequently the case instead that leading order (LO) and next-to-leading 
order (NLO) provide equally good descriptions of the data; and the NLO 
calculation is preferred mainly because it reduces theoretical systematic un- 
certainties in quantities such as ag (Q^)- As an inclusive scattering process 
at space-like momentum transfer, DIS is perhaps the phenomenon for which 
QCD is most rigorously applicable. NLO cross sections are already fully cal- 
culated and NNLO computations for several processes exist as well. Testing 
the entire theory at NNLO seems feasible. 

• Precise measurements of ag- at moderate Q"^ The running of the strong 
coupling constant is largely determined by excellent measurements of the 
r lepton hadronic decay width at s = as well as a series of precise 
measurements at s = M|. Getting this running right is important; for 
example, one of the few experimental pieces of evidence for supersymmetry 
is its success in getting the strong, weak, and electromagnetic couplings to 
unify at one scale. Given this importance, further precise measurements of 
as at scales between and m| are valuable. A z/MC provides at least two 
ways of achieving these: through scaling violations in non-singlet structure 
functions and through evaluations of the Gross-Lleweyllyn-Smith sum rule. 

• Studies of two-scale QCD via v and v heavy flavor production. A quark q is 
treated as a heavy object in DIS if Q'^ < whereas it is instead considered 

to be a parton when log (^Q^/m^^ ^ 1. Deep inelastic scattering at a z/MC 
allows study of the transition of q from heavy quark to parton by opening 
the possibility of measuring quark mass effects at a series of scales in CC 
and NC scattering. The possibihty of using vertexing targets maintains the 
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inclusive nature of measurements by avoiding the need for final state lepton 
tagging. 

• Neutrino Spin Physics. Charged lepton scattering experiments from polar- 
ized targets show that u and d type quarks carry very little of the nucleon 
spin and have hinted at strong polarization effects in gluons and strange 
quarks. A vMC creates the first possibility of using polarized targets for 
neutrinos and brings all of their power for flavor and helicity selection to 
bear on nucleon spin physics. 

• Neutrino Nuclear Physics. Thin nuclear targets at a uMC can rapidly 
acquire the statistics to make measurements of the A dependence of the F2 
structure function for neutrinos that complement those from charged lepton 
scattering. The first precise measurements of the A dependence of xF^ will 
become available as well. 



2.1 Background on Measuring Parton Distribution Functions and QCD with 
Non-Polarized Targets 



Invariance principles dictate the general form of I'M [N = p or n) nucleon 
scattering. For energies much greater than the final state lepton mass and to 
leading order in electroweak couplings: 



" ^CCjNC _ LipMjsstj^ 

dxdy ~ TT (1 + Q'lMlf 



(l-y- U O' 




(10) 



with Gp the Fermi coupling constant, Mjv the nucleon mass, E^, the neu- 
trino energy, y the inelasticity, x the Bjorken scaling variable, and the 
negative squared four- momentum transfer to the nucleon target. The plus 

(minus) sign in the final term is conventional for neutrino (antineutrino) scat- 
tering, and My = Mw {Mz) for CC(NC) scattering. The structure functions 

2^-^i!cc/ivc(^><9^), F^^c/Nci^^Q^) and xF^^^^^^{x,Q^) contain aU the in- 
formation about the internal structure of the target. The cross sections for 
electron-neutrinos and muon-neutrinos are nearly identical, up to electroweak 
radiative corrections. 

The SF depend on the A and Z of the target nucleus, on whether the beam 
is neutrino or antineutrino, and on whether the scattering is CC or NC. They 
can be experimentally extracted in principle by measuring the differential 
cross sections in fixed x and bins as a function of y and then exploiting the 
y dependences shown in Eq. 10 to fit for 2xFi, F2, and xF^. In practice, the 
reduced y coverage created by cuts on final state lepton energies and kinematic 
constraints limits this procedure, and various model-dependent alternatives 
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have been assumed. For example, the SF 2xFi has rarely been measured in 
neutrino scattering; instead it has been related to F2 through a model for the 
longitudinal structure function 



F2(x,Q^)(l + 4AfV/Q^) 
2xFi {x, g2) 



-1, 



(11) 



with Rl {x, Q'^) computed from QCD or taken from charged lepton scattering. 
Charged current interactions, with their observable lepton in the final state, 
can be much better reconstructed than NC interactions and so we will assume 
CC SF in the discussion that follows. 

Neutrino-nucleon scattering is the only DIS process that can provide mea- 
surements of the parity-violating F3 structure functions, apart from the much 
less precise measurements in a different kinematic regime from HERA. The 
parity-conserving 2xFi and F2 structure functions for neutrino-nucleon scat- 
tering probe different combinations of quarks to the analogous SFs defined for 
charged lepton DIS experiments. 

A rough summary of the current knowledge of neutrino SF follows. A more 
complete review may be found in Ref. 12. 

(1) Measurements at an accuracy of a few percent exist for n/p, v/D averaged 



using iron targets for 10"^ ^ a; < 0.7 and < 200 GeV^. The x and 
ranges are highly correlated by the limited range of beam energies. 
These measurements assume a model for Rl{x^Q'^). Uncertainties on 
F2,cc {x^ Q"^) are dominated by systematic effects, while xF^^c {x, Q"^) 
errors still contain a significant statistical contribution. 
(2) Measurements at the ~ 10% level exist for: 

(a) the n/p, vjv averaged (:X', Q^) using an iron target; 

(b) the njp averaged xF^^jq (x, Q^) — xF^^jf-, (x, Q^); 



(c) F^lc (^> Q") , Flee i^, Q') and xF^fcc, xF^Sc for k^u,P and with 



D = deuterium. 

(3) No SF-oriented neutrino experiments are currently in operation and no 
new experiments are planned other than a possible i/MC program. 

The SF goal for an z/MC is simple: to measure, over as wide a range of x and 
as possible, the six SF of Eq. 10, particularly for the proton and deuteron 
but also for other nuclear targets. 



CC SF 



F.,cc{x,Q^)=-,J: e 



F2%{x,Q^) 




k=v,v N=n,p 



(12) 
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2.2 Measurement of Quark Parton Distribution Functions 



Structure functions provide much of the information used to deduce PDFs. In 
turn, the PDFs are crucial for all predictions of event rates at the Tevatron 
and LHC. To avoid nuclear complications, z/MC SF should be extracted with 
proton and deuterium targets. 

To leading order (LO) in QCD, the SF can be expressed in terms of nucleon 
PDF as Q: 



2,CC' 



xFr{x,Q')=2 



(13) 



and X can be identified as the target's fractional 4-momentum carried by the 
struck quark. The Callan-Gross relation holds. 



vN{uN) 



^uN{vN) 



(14) 



implying that F2 
information. 



vN{uN) 



X, Q ) and 2xF^ 



vN{uN) 



X, Q"^) provide redundant parton 



Measurements of the eight independent observables of Eq. 13 (two SF for 
each beam on two targets) represents more information than is available from 
charged lepton scattering but is not enough to specify the eighteen indepen- 
dent PDF for each target N (m, d, s, c, plus their antiquarks and the gluon 



^ We adopt the convention that the PDF are given as parton probability functions 
multipUed by x, i.e. u (x, Q'^) is x times the probabiUty of finding a u quark with 
momentum fraction x. This is close to what is actually measured in experiments 
and is also the form provided in compilations such as PDFLIB[16]. 
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for n and p) . Further constraints emerge from isospin symmetry: 



7/" 


[x,Q') 


= (F{x,Q'') = 


d{x, Q^), 




{x,Q') 


= uP{x,Q') = 


u[x,Q j , 




[x,Q') 


= dP{^x,Q^) = 


d(x, Q^) , 




[x,Q') 


= uP{x,Q^) = 


u{x, Q^), 




[x,Q') 


= (x, Q') E 


E s , 




[x, Q') 


= {x, Q') E 


E s (x, Q^) , 


c" 


[x. Q') 








[x.Q') 




c(a;,Q^), 


9^ 


[x,Q') 


^9'{x,Q') = 





The nine extra constraints of Eqs. 15 reduce the number of independent PDF 
to nine. More reduction occurs if one assumes 

s{x,Q^) = s{x,Q^) , (16) 

which is supported by measurements of dimuon production in uFe and uFe 
scattering[17,18]. Furthermore, at the accessible to lower energy uMCs, 
one can consistently adopt a three-flavor QCD scheme, whence 

c (x, Q^) = c (x, Q^) = 0; (17) 

these relations are supported by measurements of NC charm production in 

Pfj^Fe and z/^-Fe scattering[19,20]. This finally leaves six independent PDF: 
u {x, (5^), d{x, Q"^), s {x, Q^), g {x, Q"^), and the valence distributions 

Uv (x, Q^) = u [x, g2) _ u [x, g2) , 

dv (x, = d (x, Q^) - d (x, Q^) . (18) 



At LO, the four deuterium SF can be combined to yield the total quark and 
valence quark distributions 



2,CC 



:x,Q') = Fl^cix,Q') 

= u{x,Q^) + d{x,Q'^) 

+ u{x, Q^) + d{x, Q^) + 2s (x, Q^) , 



(19) 



The four individual proton SF then allow separation of the u and d quark 
contributions. Sensitivity to s {x, Q"^) emerges from 



xF^gcix. Q^) - xFl^cix, Q') = 4s {x, Q') , 



(20) 
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and, in principle, from comparison to the charged lepton F2^q{x, Q'^), 



(21) 



Practical implementation of Eqs. 20 and 21 has been stymied by difficulties in 
controlling systematic errors, and current measurements of s (x, Q"^) all come 
from semi-inclusive u^N and V^N charm production. 

The gluon PDF does not enter directly at LO for CC SF. It affects the QCD 
evolution of F2 and 2xFi, appears with next-to-leading order (NLO) cross 
section terms, and enters directly into semi-inclusive double heavy quark pro- 
duction through the NC and CC processes 



as well as neutrino J/ if) and T production[21]. 
2.3 Tests of Perturbative QCD 

At a i/MC, data will be of sufficient precision to probe QCD to NLO, and 
perhaps to NNLO. 

In practice, the data will be in the form of differential cross sections in x and 
y at different neutrino energies. The general procedure for a QCD analysis of 
the data consists of: (1) defining a favorable kinematic region where pQCD is 
expected to apply to within small corrections; (2) choosing a parameterization 
of the six PDF p(x, Q^; {A}), where p — (uy, dy, u, d, s, g^, Ql is a reference 
and {A} are a set of parameters describing the PDF; (3) using pQCD to 
fit the data over all Q'^ by varying {A}. 

The quality of such a fit to the data constitutes an immediate test of QCD. 
For example, QCD can be verified to NLO if NLO pQCD provides a better 
fit than LO. While contemporary practice dictates fitting with NLO pQCD, 
it is worth pointing out that neutrino data do not convincingly favor the 
higher order calculation over a LO interpretation. Assuming a good pQCD 
fit, one can then go on to extract the single parameter of QCD itself, which 
can taken to be the strong couphng constant evaluated at the reference Q"^ of 
the experiment ag (Qq)- 

As a technical comment, the fit procedure is admittedly somewhat complex, 
and the need to treat PDF parameterizations and QCD together arises from 



(22) 
(23) 
(24) 
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kinematic acceptance issues. The value will anyway be limited to < 
2ME and cuts on final state lepton energies may further limit the available 
{x,Q^,y) phase space, particularly at lower energy i/MCs. More precisely, 
if P™™ and E^i^^^n represent the minimum acceptable final state lepton and 
hadron energies, then for a given Q^, one has 



gV2M {e - Pr) <x< min [l, /2ME^Xd 

EWXd/E <y<l- Pr/E. 



(25) 



This imphes that it is impossible to span all of x for a fixed Q'^ in order to 
extract the PDF of Eqs. 13. An interpolation scheme is needed to connect 
different regions of x and space; fortunately, pQCD provides just that 
scheme through the DGLAP equations[22-25]. 

More important even than the high statistics and potentially improved sys- 
tematics promised by i/MCs are the richness of evolution tests created by the 
availability of 12 proton and deuterium SF, even with the limited phase space. 
Examples include: 



(1) 



i [xF^cci^^Q'^) + ^Fl^ci^^Q"^)]- the "classic" non-singlet SF's evolu- 
tion is independent of g {x, Q^) but suffers uncertainties from charm pro- 
duction. 



(2) I 



x,Q' 



pvD , 



: the most precisely measurable SF use- 
ly constrains g {x, Q"^) and cross-checks charged lepton scattering. 
(3) -p2^cc(^' Q^) ~-^2^cc(^' Q^)'- ^ ^^'^ combination for a z/MC; this difference 
is both independent of g (x, Q^) and charm production. 

Sum rule-tests comprise some of the most accurately calculated observables in 
QCD. For example, the Gross Llewellyn Smith (GLS) sum rule [26,28] yields: 



GLS 



^ dx 

X 



dx 



(26) 



with a{Nf) and b{Nf) are known functions of and the specified number of 
active fiavors Nf used in the pQCD analysis. Corrections from order (as/Tr)'^ 
pQCD[27] and order M'^/Q'^ higher twist effects have been calculated[29,30]. 
A uMC would further provide the targets and the statistics to evaluate for the 
first time with a precision sufficient to test QCD through the Adler sum[31] 



(«0 



1 dx 



X 
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assuming that 2F2^; 



Structure function analyses in neutrino scattering are a complex business that 
requires a painstaking attention to systematic error sources. The current best 
QCD measurements - the evolution [32] of xF^ and the GLS sum rule [33] - are 
limited by energy calibration uncertainties in the former and flux-related er- 
rors in the latter. Calibration uncertainties will be reduced at uMCs through 
use of lower mass particle spectrometers that allow better resolution and the 
possibility of using in situ sources such as Kg and J/ip decay. The ability to 
simultaneously measure and scattering provides built in cross-checks, 
and the z/g CC scattered electrons can be measured both magnetically and 
calorimetrically to cross-calibrate the spectrometer and calorimeter. Flux er- 
rors will be diminished considerably by the simplicity of the neutrino source 
compared toir/K decay beams, which permits much more reliable monitoring. 
High rates will also allow direct flux measurement in the i^^Pg mode through 
use of the electroweak reactions i/i^e" /^"^'e and PgC^ A*~P^- 

2.4 Heavy Quark Production 

The simple language of the quark parton model must be modified for i/N DIS 
events with a charm or beauty quark in the final state in order to take into 
account the non-negligible quark mass. This presents both a challenge and an 
opportunity to test the QCD formalism for making these corrections. 

Perhaps the simplest and most widely used correction scheme is the essentially 
kinematic "slow rescaling" model of Georgi and Pohtzer [34], which amounts 
to a redefinition of the scaling variable x through 



where tuq is a heavy quark mass. One of the shortcomings of this prescription 
is that it fails to address ambiguities in the PDF that arise when one attempts 
a self-consistent pQCD treatment of heavy quark production. 

Neutral current charm production illustrates this ambiguity. On the one hand, 
this can be treated as a "flavor-creation" process, Z + g cc, which occurs 
at order a] in the pQCD expansion. In this case, the appropriate scheme for 
calculation of QCD radiative corrections involves only three light flavors (the 
so-called three-flavor scheme) , and all of the effects associated with the charm 
quark arc accounted for perturbatively in the hard scattering. On the other 
hand, the charm quark can be considered to be "light" at high energies when 
3> so the mass of the charm quark is much smaller than the relevant 
physical scale. In this case, one can view the charm quark as a parton with a 




(27) 
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corresponding parton distribution function (the so-called four-flavor scheme). 
In particular, the lowest order Wilson coeflicient for charm production in this 
situation is a "flavor-excitation" process, Z + c —>■ c, which is of order a^. 
These two viewpoints must be reconciled in a consistent manner in order to 
avoid double counting. 

The corresponding calculational problem for pQCD is the presence of a heavy 
quark mass scale that can be comparable to the of the interaction. In the 
absence of this complication, the factorization theorem of QCD separates the 
high scale set by the value of transferred momentum, /i^ ~ Q^, from the 
low scale of hadronic physics, /i| ~ ^qcd^ operator product expansion 
techniques can be used to sum the large logarithms of the form \og[ii\/ jil) 
that multiply the expansion parameter, ag- 

A practical recipe for incorporating heavy quark masses in the pQCD summa- 
tions is provided by the ACOT [35] prescription, which treats the number of 
active flavors, A^/, as a scale-dependent quantity. In particular, it suggests the 
use of a three-flavor evolution for fi < rUc and a four-flavor evolution above 
mc, with continuity at the break point. In this prescription, the parton dis- 
tribution functions are labeled by the number of active flavors and the heavy 
quark parton distribution functions, /Q/Ar(C,/x), vanish for [Xh < rnq and sat- 
isfy the usual MS QCD evolution equation (with massless kernel functions) 
for iih> mq. 

Charm and beauty production at vMCs will include large numbers of events 
with above, below and around the effective scale and should provide the 
most stringent tests from any experimental process of the pQCD formalism for 
heavy quark production. At that point, the experimental precision will likely 
require a more careful quantum field theoretical treatment of heavy quark 
masses. 

Available rates at a z/MC (Figs. 4 and 5) provide access to five different values 
for itlq: rric (from UfN — > frcX), 2mc (from UfN UiCcX), mi, (from UfN 
i~bX), rrih + rric (from u^N — £^d)X), and 2mh (from u^N — UfbbX). To 
the degree a consistent pair of values for m?, and rric emerges, one would 
have established theoretical control over these production processes. Those 
mass values could also be compared to what QCD sum rules, heavy quark 
expansions and lattice QCD yield. They can also be used to calibrate models 
for heavy quark production that enter into other z/MC analyses, such as the 
pQCD tests described above and precision sin^ 6w measurements in inclusive 
NC scattering discussed in Sec. 4. 

Other than the ubiquitous high rates, the main experimental asset for a z/MC 
in these studies will be the opportunity to use vertexing targets. This should 
allow for lifetime tagging of heavy flavors that minimizes systematic effects 
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due to fragmentation and decay uncertainties. 

2. 5 Parton Distribution Functions at Large Bjorken x 

The unique level of quark-by-quark characterization of nucleon structure ex- 
pected at i/MCs will provide an invaluable reference source for many diverse 
analyses in collider and fixed target physics including, of course, other precision 
analyses at z/MCs. Precise measurements as a; — > 1 are particularly relevant 
to the modeling of rates for interesting physics processes and backgrounds at 
hadron colliders because uncertainties at high x and at the typical val- 
ues for i/MCs will evolve to uncertainties at much lower x as increases to 
collider values. Uncertainties at high x in current nucleon PDF derive from 
two sources: the ratio d{x, Q'^)/u{x, Q'^) as x — > 1 and the role of higher twist 
corrections. 

Analyses on present leptoproduction data sets that used hydrogen and deu- 
terium targets have been unable to precisely pin down the high-x behavior of 
d{x)/u{x). For example, QCD fits [36] to the high-x NMC data and the CDF 
W^-decay asymmetry improve with the inclusion of a simple x-dependent cor- 
rection for a d{x)/u{x) ratio that asymptotically approaches 0.2 rather than 
zero for X — > 1. On the other hand, the recent CTEQ5 global QCD analysis 
[37] found that including this correction had httle effect on the quality of their 
fits. 

Besides the statistical and experimental uncertainties in the data sets, a com- 
plication with high X analyses is the need to model nuclear binding effects in 
deuterium This issue can be avoided at a z/MC, where a high statistics expo- 
sure to a H2 target alone could directly measure the d{x)/u{x) ratio in protons 
as a; — > 1 from the ratio of to u^p cross sections. Such a measurement would 
require only a small correction for the residual sea quark contributions at high 

X. 

Measurement of quark PDFs at high x is closely related to the question of the 
leading power corrections in the QCD perturbative expansion that are known 
as "higher twist effects" . The n*^ order higher twist effects are proportional to 
l/Q^" and reflect the fact that quarks have transverse momentum within the 
nucleon and that the probe becomes larger as decreases, thus increasing 
the probability of multi-quark participation in an interaction. As was the case 
with the u/d ratio, different analyses of higher twist corrections in current data 
leave unresolved issues that would benefit from new experimental information. 

An analysis by Milsztajn[38] that combined electroproduction data from SLAC 
with BCDMS muo-production data found that the relative size of the twist- 
4 contribution rapidly increased above x — 0.4 and was equal in magnitude 
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to the leading l/logQ^ term for x > 0.75. The only measurements of this 
higher-twist term in neutrino experiments have been two low-statistics bubble 
chamber experiments: in Gargamelle [39] with freon and in BEBC with NeH2- 
Both bubble chamber analyses are comphcated by nuclear corrections at high- 
X. However, both found a twist-4 contribution that is smaller in magnitude 
and, most significantly, of opposite sign to that of the charged leptoproduction 
analysis. 

In contrast, a CTEQ global QCD analysis that combines neutrino and charged 
lepton production analyses up to a; = 0.75, finds that no higher-twist term is 
required for values of Q'^ down to 0.7 GeV^. However, this analysis uses a cut 
on the invariant mass of the hadronic system, W > A GeV, that could exclude 
the bulk of any twist-4 contribution. Most recently, the Yang-Bodek analysis 
mentioned above [36] reanalyzed electroproduction data looking for a higher 
twist contribution. They find that by incorporating a NLO QCD analysis, as 
opposed to the LO analysis used by Milsztajn, the higher twist contribution 
becomes much smaller. However, they have also included the d/u — > 0.2 model 
in their extraction of higher twist, which will reduce the size of any extracted 
higher twist term. 

From a more theoretical viewpoint, a recent CTEQ paper by Guo and Qiu [40] 
has predicted the leading x dependence to the higher- twist term in F2(a;, Q"^) 
in the region of large x and found that the higher- twist contribution is different 
for u and d quarks. Were this to be true, then the d{x)/u{x) analysis and the 
leading power correction analysis would be directly intertwined. Interestingly 
enough, this prediction only depends on two non-perturbative parameters and 
is therefore highly constrained. It could be tested at both CEBAF and a uMC 



2. 6 Examining the Spin Structure of the Nucleon 



A unique new feature of vMCs would be the availability of sufficiently intense 
beams to allow, for the first time, measurements of neutrino scattering off nu- 
cleons in polarized targets. This would provide access to the polarized nucleon 
structure functions for neutrino- nucleon scattering and could answer several 
currently unresolved questions about the spin structure of the nucleon. 

Polarized structure functions (PSFs) can be represented approximately by the 
quark-parton model in terms of the differences between the parton densities 
of quarks polarized parallel to the nuclear spin and those that are polarized 
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anti-parallel: 



5q{x,Q^) = q^\x,Q^) - q^\x,Q^), (28) 

5q{x, Q^) = q^\x, Q^) - q^\x, Q^).5g{x, Q') = g^\x, Q') - g^^{x, Q'). 

(29) 



In the naive (LO) quark-parton model, the polarized structure functions gi 
and g^ have quark spin content corresponding to the quark content of the 
parity conserving and parity violating unpolarized structure functions Fi and 
F3, respectively: 

gf(x, g2) = sd^{x, + ss^'ix, g2) + su^'ix, q^) + Sc^'ix, q^), (30) 

and 

2xg^/{x, g2) = - [6d{x, g2) + 6s{x, Q^) - 5u{x, Q^) - 5c{x, Q^)] , (32) 
2xgf{x, g2) = -[6u{x, g2) + Sc{x, Q^) - 6d{x, Q^) - 5s{x, Q^)]. (33) 

The other parity-conserving PSF, g2., has no simple interpretation in the LO 
quark-parton model, as is also the case with the parity-violating g^^. The re- 
maining PSF, 5f5, is related to g^, via g^ = 2xg3 in the naive quark-parton 
model. 

The levels of the different polarized quark PDF's can be regarded as quanti- 
fying the extent to which a parton of flavor q "remembers" the polarization 
of its nucleon parent in interactions at a given Q^. The nucleon spin, (|), can 
be decomposed as: 

i = i (AE(g^) + Ag{Q') + L,{Q') + L,{Q')) , (34) 

where AS = Au + Ad + As + . . . and Ag are the integrated net quark and 
gluon helicities along the nucleon spin direction, e.g., 

Au(Q'^) = C 5u{x, Q^)dx, 
Jo 

Ag{Q^)= f' Sg{x,Q^)dx. (35) 
Jo 

The Lq and Lg represent the relative orbital angular momenta of the quarks 
and gluons, respectively. 

To date, the only charged lepton polarized DIS data comes from non-collider 
experiments, running at energies where photon exchange dominates. Thus, 
only the parity conserving polarized structure functions gf^ and 512'^ have 
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been studied so far. The resTilts are surprising [?]. From the measured first 
moment of g['^ structure function, it can be concluded that only 30% of the 
nucleon spin is carried by quarks if one assumes Ag = 0, i.e., AE = 0.3. 
An additional conclusion in this simple parton model interpretation is that 
the strange sea is anti-aligned with the nucleon spin, As = —0.1. In the 
more realistic QCD-enhanced parton model, the QCD evolution of the quark 
distributions brings in contributions from gluon radiation (g qg) and pair- 
production {g — > qq) that induce a non-zero gluon spin distribution. This 
latter conclusion is supported by pQCD analyses of the evolution of the 
available data on g{^ , as is illustrated in Fig. 6. 

Although the decomposition of Fig. 6 is highly model dependent, the abil- 
ity of neutrinos to select specific quark flavors should allow a polarized DIS 
experiment at a muon storage ring to greatly constrain such models. 

2. 7 Experimental Setup and Measurement Technique 

An ideal polarized target would have a mass of at least 10 kg, high polarizabil- 
ity, a large fraction / of polarizable material, and the capability for polarizing 
both protons and neutrons. The sohd HD compound material "ICE" appears 
to be a promising new technology [41]. Both the H and the D can be polarized, 
either in separate experimental runs or together. The expected polarization 
and dilution are Ph = 80% and /^^ = 1/3 for hydrogen, and Pd=50% and 
fi) = 2/3 for deuterium. A 7 kg (density pt=0.11 g/cm^) polarized target with 
the qualities mentioned above can be built out of an ICE-target that is 20 cm 
in radius and 50 cm deep, perhaps mounted upstream from a general purpose 
neutrino detector like that in Fig. 2. Raw event rates in such a polarized 
target could be of order 10^ events per year, as shown in Table 2. A previous 
study for uMCs [42] has shown that a 200 kg target of polarized butanol with 
10% polarization could measure the strange sea polarization to about 3% in a 
one year run in a beam downstream of a 250 GeV muon collider with a 10 m 
straight section. 

Measurement of gi requires both i> and V beams and follows from the double 
cross section difference: 

where a'^^ denotes, for example, the neutrino scattering off a target N with 
its spin polarized parallel to the neutrino helicity. This difference should be 
measured as a function of x and y to test for effects of unwanted PSF such as 
g2. For g^ one evaluates the similar quantity 
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Fig. 6. Model-dependent decomposition of nucleon spin into contributions from 
quarks and gluons. The two vertical dotted lines show the naive QPM expectation 
and the results from a NLO fit to most of the available data on g[. 

Cross checks can be made with the target transversely polarized, where the 
suppression factors for the unwanted PSF are different. 

2. 8 Applications of Polarized Parton Distribution Data from vMCs 

Measured values for gi from any nucleon target at a z/MC would allow verifi- 
cation of the predicted decomposition: 




(38) 
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where the factor C is model dependent and a common choice is C = Nfas/27r, 
and one has knowledge of AS from other data. In contrast to polarized charged 
lepton scattering, this determination can be done without relying on low en- 
ergy input from beta decay data augmented by SU(3) symmetry. 

For another application, one notes that g^, and (75 probe only non-singlet 
combinations of parton densities so that they do not get a contribution from 
the gluon density. Then, for scattering off an isoscalar target: 

xgl^ - xgf^ ^Sc + Sc-Ss-Ss . (39) 

(This is the polarized target analog of Eq. 20.) Assuming that Sc -C Ss, 
this provides a measure of the level of polarization of the strange quark sea. 
Like xF^, the non-singlet PSF, g^, g^ and g^, have QCD evolutions containing 
no contribution from gluons at lowest order. Comparison of the non-singlet 
functions with the singlet SF gi and F2 should therefore provide an indirect 
means for measuring the gluon contribution A^f. 

As was the case for unpolarized targets, final state quark flavor tagging should 
provide additional semi- inclusive structure functions to augment the informa- 
tion on quark content from inclusive polarized structure functions. Probably 
the most important example will be studies of the strange spin contribution 
via charm production. Figures 4 and 5 show that approximately 5% of the 
events from a 50 GeV muon storage ring will have charm in the final state, 
with 20% of these decaying semi-leptonically. An experiment with 20 million 
neutrino interactions thus will have 2 x 10^ semi-leptonic charm events before 
kinematic cuts, which should be sufficient for a precise measurement of 5s and 
5s. 

2.9 Studying Nuclear Effects with Neutrinos 

Nuclear effects in DIS have been studied extensively using muon and electron 
beams but have only been glanced at for neutrinos in low-statistics bubble 
chamber experiments. Neutrino experiments with high statisticss have for- 
merly only been possible using heavy nuclear targets such as iron calorimeters 
and, for these targets, nuclear effects in vN interactions have typically been 
considered as problems to overcome rather than as a source of physics insights. 
This section reviews the physics of nuclear QCD that is relevant to neutrino 
interactions and shows that i/MCs could instead provide experimental condi- 
tions where a great deal of interesting knowledge could be added to this field 
by using a variety of heavy nuclear targets as well as H2 and D2. 

Nuclear studies at a z/MC could use a general purpose detector such as that 
in Fig. 2 if it could be designed with interchangeable targets. Alternatively, 
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a detector dedicated to nuclear studies could be used, perhaps with a geome- 
try similar to that of the Fermilab E-665 Tevatron muon experiment [43] . For 
example, it could consist of a liquid H2 or D2 target followed by a rotating 
support of targets with different A interspersed with tracking chambers, and 
then by an appropriate calorimeter/muon spectrometer. A i/MC should be 
easily capable of supplying the event sample sizes required to examine the 
predicted nuclear effects. For example, one can consider a one-year exposure 
of each target to the beam in the 50 GeV i/MC scenario of Table 1. Then 10^ 
events would be acquired in targets subtending out to the 1/7^ characteristic 
angular size of the neutrino beam and that had lengths of 18 cm for D2, 1.4 
cm for graphite and only 0.16 cm for tungsten. 

Different types of nuclear effects arise on passing through four distinct regions 
in Bjorken x: 

(1) "shadowing" for x < 0.1, 

(2) "anti-shadowing" for 0.1 < a; < 0.2, 

(3) the "EMC Effect" for 0.2 < x < 0.7, and 

(4) "Fermi motion " for x > 0.7. 

These regions will now be discussed in turn. 
2.9.1 Low x: PC AC and Nuclear Shadowing 

In the shadowing region, x < 0.1, there are several effects where a neutrino 
probe could provide different insights to charged lepton probes. Considering 
first the limit as 0, the vector current is conserved and goes to zero but 

the axial-vector part of the weak current is only partially conserved (PCAC) 
and F2{x,Q'^) approaches a non-zero constant value as — > 0. According 
to the Adler theorem [31], cTj^tv can be related to to a^^N at = 0. A i/MC 
should be able to address the question of what effect a nuclear environment 
has on the Adler theorem. 

The region of vector meson dominance (VMD) is reached in nuclear scat- 
tering of charged leptons {i'^A scattering) in the \ow-x shadowing regime as 
increases from but remains below of order 10 GeV^. The physics con- 
cept of VMD is the dissociation of a virtual photon into a q — q pair, one 
of which interacts strongly with the "surface" nuclcons of the target nucleus. 
Thus the interior nucleons are shadowed by the surface nucleons. Neutrino 
scattering should involve not only a VMD effect (though now with dissocia- 
tion of a virtual W rather than a photon) but also additional contributions 
from axial- vector mesons such as the Ai. Further, there should be additional 
non-perturbative effects that appear as nuclear shadowing (mainly in large 
nuclei) and which involve gluon recombination from nucleons neighboring the 
struck nucleon that shift the parton distributions towards higher values of x. 
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Boros et al. [44] predict that the resulting shadowing effects in V(,A scattering 
will be roughly 1/2 that measured in i^A scattering. A more quantitative 
analysis by Kulagin[45] uses a non-perturbative parton model to predict that 
shadowing in V(_A scattering to be either equal to or slightly above that in 
t^A scattering at = 5 GeV^. It also attempts to determine the quark 
flavor dependence of shadowing effects by separately predicting shadowing for 

F^^cM^). F^NcM) and xF^^cM) ■ 

The predictions of Kulagin should be testable at a i/MC, as is indicated by 
Fig. 7. This shows a simulation for a z/MC data sample of the predicted mea- 
sured ratios of the values for F2{x) and -^3(2^) on a calcium target divided by 
those on a deuterium target: 

R2a{x) = F^ix, Q')[Ca]/F2ix, Q'')[D2], (40) 
Rsa(x) = Fs(x,Q^)[Ca]/Fs(x,Q^)[D2], (41) 

where the ratios are an event- weighted average over the experimental range. 
As can be seen, the predicted difference between the shadowing on sea and 
valence quarks is clearly visible down to x — 0.02 — 0.03. 

2.9.2 Mid x: Anti- shadowing and the EMC Effect 

Drell-Yan experiments have also measured nuclear effects and their results are 
quite similar to DIS experiments in the shadowing region. However, in the 
anti-shadowing region where R2A niakes a brief but statistically significant 
excursion above 1.0 in DIS, Drell-Yan experiments see no effect. This could be 
an indication of difference in nuclear effects between valence and sea quarks. 

Eskola et al. [46] have quantified this difference using a model which predicts 
that the differences between nuclear effects in xF^QQipCiQ"^) and Fl^^fjipCiQ"^) 
should persist through the anti-shadowing region as well. Taking the work 
of Kulagin and Eskola together implies that nuclear effects in xFi^^Q^x.^Q'^) 
should be quite dramatic, with more shadowing than F2cc{x,Q'^) at lower x 
followed by i?3^ rising fairly rapidly to yield significant anti-shadowing around 
X = 0.1. The 50 GeV uMC experiment assumed for Figure 7 should be able to 
measure anti-shadowing effects and the difference between shadowing effects in 
F2{x, Q"^) and xF2,{x, Q"^) at a statistical level of about 6 standard deviations. 

2.9.3 High x: Multi- quark Cluster Effects 

Analyses from DIS experiments of F2{x,Q'^) in the "Fermi-motion" region, 
X > 0.7, have required few-nucleon correlation models and multi-quark cluster 
models to fit the data. These models boost the momentum of some of the 
quarks, producing a high-x tail in F2{x,Q'^) that is predicted to behave as 
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Fig. 7. Simulated predictions for the ratios of both the F2 and xF^ structure 
functions between calcium and deuterium, plotted as a function of Bjorken x and 
assuming the theoretical model of Kulagin [45] and data samples from a z^MC. The 
sizes of the error bars correspond approximately to 1-year exposures for both targets 
at the 50 GeV uMC of Table 1. 



e However, fits to fiCa [47] and uFe [48] have obtained two different 
values for the fitted constant a: a^ca = 16.5 ± 0.5 and a^pe = 8.3 ± 0.7 ± 
0.7, respectively. This was considered surprising because of the expectation 
that any few-nucleon-correlation or multi-quark effects would have already 
saturated by carbon. A high statistics data sample from a z/MC, with precise 
kinematic reconstruction of the Bjorken x values for each event, could go a 
long way towards resolving the dependence of the value of a on the nucleus 
and on the lepton probe. 



3 Studies of the CKM Quark Mixing Matrix 
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3. 1 Introduction 



The CKM quark mixing matrix comprises a set of fundamental parameters 
of the Standard Model that reflect the dynamics shaping the generation of 
particle masses. Furthermore, it provides essential input for predictions on 
various CP asymmetries in B decays which should be sharpened as much as 
possible in order to increase their sensitivity to the presence of new physics 
processes. 

It is likely that i/MCs could make a very substantial contribution to measure- 
ments of quark mixing, with systematically unique measurements of four of 
the nine moduli in the CKM matrix that should be comparable to or better 
than the best future measurements at other experiments. The relevant moduli 
are \ Vcd\-i |V^s|, and |V^b|, to be extracted from charm production off d and 
s quarks, and from beauty production off u and c quarks, respectively. If muon 
colliders ever reach the 100 TeV energy center-of-mass range then their z/MCs 
could further provide a unique opportunity for precise direct measurements of 
the CKM elements involving the top quark [14] . 

This section is organized as follows. Section 3.1.1 gives a summary of the 
current experimental status and future expectations for before the advent of 
i/MCs . Section 3.1.2 then gives a general overview of the method for extract- 
ing CKM information at i/MCs . Sections 3.2 and 3.3 then give more specific 
details on vMC experimental analyses involving the production of charm and 
bottom quarks, respectively. A further analysis that may have potential for a 
measurement of CKM parameters - diffractive charmed vector meson produc- 
tion - is briefly addressed in subsection 3.4. Finally, subsection 3.5 summarizes 
all of these expected improvements from i/MCs to our knowledge of the CKM 
matrix. 



3.1.1 Current Experimental Knowledge of the Relevant CKM Matrix Ele- 
ments, and Future Expectations 

Table 3 presented the current [15] experimental values and uncertainties for 

the absolute squares of the interesting CKM matrix elements for z/MCs : Vcd, 
Ksj Kib and Vc5. The first two of these can be more tightly constrained by 
adding the assumption of three family unitarity, so improved measurements 
of these elements can be considered as testing these unitarity conditions. 

The values for \Vcd\ and \Vcs\ have been extracted from charm production in 
deep inelastic scattering supplemented by information gathered in semilep- 
tonic D decays, while |V"cb| and \Vuh\ have been obtained from semileptonic B 
decays. 
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In the next few years, exclusive semileptonic charm decays hke D — > IvKj 
K* j n / p wiU be measured more precisely. However it is not clear at present 
whether theoretical technologies like QCD simulations on the lattice will im- 
prove sufficiently to fully utilize such data for extracting more accurate values 
of \Vcs\ and \Vcd\- 

At present |Vc6| is measured with about a 5% theoretical uncertainty and a 5% 
or better experimental uncertainty. (For comparisons with the uncertainties 
quoted in Table 3 recall that the fractional uncertainty in the square of the 
modulus is twice that of the modulus itself.) The two most reliable determi- 
nations are from the semileptonic B width and the form factor for B iuD* 
at zero recoil. One can reasonably expect the theoretical uncertainty in these 
processes to go down to about 2% [50]. 

The situation is considerably less satisfactory for \Vub\- Most of the analyses 
employed so far contain a large dependence on theoretical models, the accuracy 
of which is hard to evaluate. One can hope that studies of both exclusive and 
inclusive semileptonic B decays done at the B factories will yield values for 
\Vub\ with a theoretical accuracy of no more than 10% by 2005. However, this 
is not guaranteed. 

It is highly desirable to improve on this situation by obtaining more precise 
and reliable values for these CKM parameters. The main reason is that they 
are theoretically expected to be a consequence of the generation of quark 
masses which in turn reflects dynamics operating at presumably ultra-high 
energy scales. It turns out that qualitatively quite distinct "textures" for the 
Yukawa couplings assumed to apply at GUT scales lead to CKM parameters 
that, due to differing renormalization effects, are numerically relatively similar 
at electroweak scales. Secondly, various CP asymmetries in B decays will be 
measured presumably to better than 5% over the next ten years; and their 
predicted values depend crucially on the size of Vub, Vch (and Vtd)- To exploit 
the discovery potential for new physics to the fullest one wants to match up 
experimental and theoretical accuracy. 

To be more specific: 

• One wants at least to confirm the values for | V^^j and \ Vch\ by a systematically 
different method. 

• One would like to determine a precise value for \ Vcd\ without imposing three- 
family- unitarity. 

• It is an important goal to extract the value of \Vuh\ with considerably less 
than 10% uncertainty. 

It appears that z/MCs are up to these tasks and possibly more. Indeed, if 
a future i/MC in the 100 TeV energy range were ever to be realized, then 
even top quark production could be studied and this would uniquely allow the 
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extraction [14] of Vta and Vts in a relatively straightforward way! 

3.1.2 Extracting CKM Matrix Elements in vMCs : an Overview 

The experimental quantities used for all the quark mixing measurements in 
z/MCs are the production cross sections and kinematic distributions of heavy 
flavor final states from CC interactions. These are related back to the quark 
couplings through a parameterization of the initial state quark distributions 
and a model of the production processes that includes threshold suppressions 
due to the quark masses. Measurements of CKM matrix elements at i/MCs 
will be analogous to, but vastly superior to, current neutrino measurements 
of iKdp that use dimuon events for final state tagging of charm quarks, which 
are reviewed elsewhere[12]. This improvement is due to the overall higher rate 
of interactions, the ability to tag charm inclusively by vertexing, and the 
capability of fully reconstructing the event kinematics for charm decays to 
hadrons. 

In principle, the struck quark can be converted into any of the three final 
state quarks that differ by one unit of charge. In practice, production of the 
heavy top quark is kincmatically forbidden except from neutrinos at energies 
above about 16 TcV, and the production of other quark flavors is influenced 
by their mass. After correcting for this kinematic suppression, the Standard 
Model predicts the probability for the interaction to be proportional to the 
absolute square of the appropriate element in the CKM matrix. The most rele- 
vant theoretical uncertainty lies in the treatment of the production thresholds. 
Measuring these cross sections at different energies and separately for CC and 
NC reactions should help in understanding the threshold behavior. 



3.2 Analyses Involving Charm Production: the Extraction ofVcd and Vcs 

The extraction of \Vcd\ will be the cleanest of the four CKM measurements at 
uMCs discussed here and is the one CKM element whose modulus is already 
best determined in neutrino-nucleon scattering. It is measured from charm 
production off valence d quarks in isoscalar targets, and this valence quark 
distribution can be accurately determined from nucleon structure function 
measurements. 

Figures 4 and 5 of Sec. 1 show that charm production occurs in ~ 7% of the CC 

interactions from 100 GcV neutrinos, and will occur at slightly lower/higher 
rates for lower/higher energies according to the changing level of mass thresh- 
old suppression. This implies that realistic charmed event sample sizes might 
reach the 10^ level. 
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The feature of the \Vcd\ analysis that makes it systematically rather clean is 
that charm production from valence quarks occurs for neutrinos but not for 
antineutrinos. Valence quarks have a harder x distribution than sea contri- 
butions from both s and (at a lower level) d quarks. It is also helpful that 
the s and d seas seen by neutrinos are also closely equal to their antiquark 
counterparts that are probed by antineutrinos, as would be verified in struc- 
ture function analyses of the uMC data set ( Sec. 2). This allows use of the 
u sample as an effective background subtraction for the sea contributions to 
the charm production by neutrinos. 

To illustrate, one can measure scattering from an isoscalar target D with u/V 
beams and extract the ratio 



dcr'^^x I dxdy — da'^^x I dxdy 
da^'^ldxdy - da^^jdxdy 



E) = .!^n;..,. . .Jn (42) 



where da^J^x^'' /dxdy represent the charm production cross sections, and (ic^^c'^^V dxdy 
denote the total inclusive CC cross sections. To leading order in QCD, this 
ratio can be computed from 



1 _ {q[x,Q')-q{x,Q^)){\V^,\'-{l-yf) 



T^o (a;, E) \VaV {q (e, Q') - q (C, Q')) (1 - 1/ + xy/0 



+ 1, (43) 



where q{x,Q'^) = {u{x,Q'^) + d{x,Q'^)) /2, ^ = x {1 + ml/Q'^), and rric is 
the charm quark mass. Note that r^^ {x, y, E) depends only on the precisely 
measured \Vudf , the well understood valence quark distribution, q{^,Q^) — 
q (^) Q^) ) the charm mass, and {Vcdl^. Because the valence distribution peaks 
at relatively high x, threshold effects associated with rric are minimized. Next- 
to-leading-order calculations are more complicated but leave these features 
intact. 

As in analyses of today's neutrino experiments, the d and s contributions to 
the charm event sample will in practice be separated from one another in a 
fit involving the measured x distributions of the events and with the charm 
quark mass as a fitted parameter. With such huge statistics, the statistical 
uncertainty in the measurement will be only at the 10~^ level so the measure- 
ment will be dominated by systematic uncertainties. Major contributions to 
the uncertainty are likely to be: 

• Estimation of the charm tagging efficiency. This error is minimized by the 
clean, high efficiency charm tagging that is possible with high performance 
vertexing and a negligible B hadron background [9]. 

• The charm production model. It is a very attractive feature that the mea- 
surement is effectively of the ratio of charm production to the total produc- 
tion in each kinematical bin. Therefore, one is sensitive only to the kinemat- 
ical suppression of charm states and not to the initial parton distribution 
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functions. Mass threshold corrections that are apphed according to theo- 
retical models can be checked and calibrated from the trends in the data 
itself over the range of kinematical bins. Therefore, this uncertainty should 
initially decrease with increasing sample size. 
• Estimation of the extent to which the sea is symmetric in quarks versus 
antiquarks. This uncertainty is lessened because the signal is largely at 
higher x than the background. If helpful, a cut in x value could be applied 
to the fit to further reduce this uncertainty. 

With experimental handles on the major systematic uncertainties, it may be 
guessed that the irreducible theoretical uncertainty on IK^P due to violations 
of quark-hadron duality might be at the percent level. 

A measurement of iKsp presents more challenges because this matrix element 
always appears in combination with the strange quark PDF s {x, Q"^), which is 
difficult, but not impossible, to measure separately. Perhaps the best scheme^ 
involves a simultaneous analysis of inclusive charged current scattering and 
semi-inclusive charm production [51] on deuterium targets. The inclusive CC 
cross section is of the form 



, uD(uD) 

dxdy 



1 yD(vD) 

dxdy 



+ 



da 



pD(vD) 

dxdy 



(44) 



where da ^^pc/ dxdy represents the part of the cross section with no charm in 
the final state {h production is neglected). To leading order 



TT 



2G'pME dxdy 

TT 



2Gj,ME dxdy 

TT 



\Vcsfs{^,Q') + \Vj'q{^,Q' 
\Vcsfs{^,Q^) + \Vj'q{^,Q' 



2GlMEdxdy .^-l' ^ ^ V l^"^"' K^' 
+ g(x,g2) {l-y)\ 



mz 



m 



2 \ 



2MEi) 



(45) 
(46) 

(47) 
(48) 



Measuring da^^x I dxdy and da'^^^x/dxdy gives da^^^/dxdy. The y dependen- 
cies of the non-charm cross sections allow measurement of q (x, Q^) and q (x, Q'^) 
Assuming that {Vudt and \Vus\^ contribute negligible error, and that s (x, Q"^) = 
s (x, Q"^), one then extracts s (x, Q"^) from the Cabibbo suppressed s — > m pieces 



of da'j^l^ / dxdy , and then \Vcs\^ (and iKdH from da^^^x / dxdy . 



Methods to extract s (x, Q^) discussed in Sec. 2 assume a known CKM matrix. 
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This procedure will require very high statistics and extreme care with theo- 
retical and experimental systematic errors as it turns on the \Vusf s {x,Q'^) 
term in the cross section, which contributes only ~ 0.15% of the total cross 
section. Inclusive and charm cross sections must be cross-normalized to high 
accuracy, and the expressions (45-48) must be generalized to include higher 
order QCD and possible non-perturbative effects. 

Additionally, it may also be possible to measure the total strange quark con- 
tent of the nucleon in NC "leading particle tagging" with a 0(1020) final 
state [9], although no experience with this method yet exists in neutrino 
physics. 

In the end, it will likely be difficult to compete with the indirect constraint 

I 1 2 h 

on \Vcs\ from precise measurements of the hadronic width of the W, , at 
LEP2, the Tevatron and the LHC. This observable is, to leading order, 

Q=u,c q=d,s,b 

where F^* is the total hadronic width. Because |T4&|^ is small and all other 
terms in the sum can be measured precisely elsewhere, a strong \Vcs\ con- 
straint emerges. 



3.3 Analyses Involving Bottom Production: the Extraction ofVub and Va, 

Measurements of Vuh and Vch from the two transitions with a h quark in the final 
state - u ^ h and c ^ b, respectively - require analyses that are conceptually 
similar to those for charm discussed in the preceding subsection. Figures 4 
and 5 show that statistics for CC b quark production to be of order 10^ events 
for a total sample of 10^*^ inclusive events, depending on the energy of the 
muon storage ring and the consequent threshold suppression due to the b 
quark mass. Besides the threshold suppression, the main reasons why the B 
production levels are so low for the two processes are that, in the one case, the 
production from u quarks is suppressed by |T46p ~ O (10~^), while production 
from charm quarks is inhibited by — 1-6 x 10~^ and the extra mass of 

the be final state. 

Given the nearly optimal vertexing geometry possible at a uMC , separating 
out much of the fractionally small b hadron event sample [9] from the charm 
background should be feasible. Assuming that the B tagging efficiency and 
purity are well known, then the b production analysis should be a relatively 
straightforward copy of that for charm, involving a simultaneous fit to the u 
and c contributions to B production from neutrino and antineutrino scattering. 
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In analogy to charm prodTiction. the resolving power of the fit arises because 
antineutrino (but not neutrino) interactions give a high x contribution from u 
valence quarks, while c ^ b transitions will be equal for u and u and typically 
at lower x. It is again easy to see the effect at LO in QCD, as follows. 

Define, in analogy to Eq. 42, 

^' - daf^/dxdy - dafjdxdy ' 

where now da^^x^"* / dxdy represents the beauty production cross sections. 
QCD then predicts 



1 {q{x,Q^)-q{x,Q^))[\V^,\^~{l-yf) 



(x, y, E) IK,!' {q (C', Q') - q (f, Q^)) (1 - I/) (1 - xy/^') 



+ 1, (51) 



with ^' = x{l-\-ml/Q'^) and the h quark mass. As was the case with 
charm, r^'~' {x, y, E) depends only on the high x valence quarks, which helps 
reduce the more substantial suppression associated with the higher value of 

Reasonably accurate estimates of the measurement precisions would require 
both a more detailed study and a knowledge of the expected event sample size. 
The precision of the |Ki5p measurement may approach the statistical limit of 
around 1% for a total neutrino event sample of order 10^° events. 

The c — > 6 transition can (unlike s — > c) be identified experimentally by 
the presence of soft c quark observed in association with the h quark in v 

scattering. At z/MC energies, it is reasonable to assume that he production 
will be dominated by 14^— gluon fusion, and a good knowledge of the gluon 
PDF at high x will be required. This may well limit the iKbp measurement 
accuracy to the few percent level. 

It is widely expected that \Vch\ will be known from B decays in other experi- 
ments to a few percent or better by 2005 [50]. Given this, an alternative and 
maybe more useful analysis strategy for a z/MC might be to combine the Vcb 
analysis data at the z/MC with the value of Vcb as extracted from semileptonic 
B decays in other experiments, with the purpose of analyzing the mass sup- 
pression in beauty production close to threshold, or in constraining the gluon 
PDF. This information could then be used as input to the extraction of Vuh 
as sketched above. 
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3-4 Vcd/Vcs via Diffractive Charmed Vector Meson Production 

It is possible to diffractively produce charmed vector mesons via VF— boson- 
pomeron scattering[52,53]: 

i^i^A {W+V) A I2-D*+A, (52) 

Ui,A (W+V) A ii-D*+A. (53) 

The Dp' cross section is of order 0.01 fb/nucleus, while that for D*^ should 
be smaller by a factor of iKd/Ksl^- For 10^° CC events, one expects of order 
1.5 X 10^ events and 7500 D* events. If the W — Dg and W — D* dynamical 
couplings are identical, the relative rate can measure the ratio of the CKM 
matrix elements. 

The experimental signature is fairly unique: a two- or three-prong muon-meson 
vertex with a / or D~g secondary. The cross section will peak at low 
momentum transfer and there should be no evidence for nuclear breakup. 

The chief theoretical uncertainty likely has to do with evaluating SU(3)- 
breaking effects in the W^— vector meson couplings and this may limit the 
potential for CKM studies from this process. Again, an alternative perspec- 
tive can be taken: to use the information available on to analyze the 
SU (3) pattern of diffractive D* vs. production. 

3.5 Improved Knowledge of the CKM Matrix from uMCs 

The Particle Data Group[15] assigns the following uncertainties to the four 
CKM parameters that were the main focus of our discussion: 



AKd|no«; = 7% [3%] (54) 

|AK..|no«, = 15% [2%] (55) 

^Vub\now = 30% (56) 

lAKblno^; = 8% (57) 



where the numbers for Vcd and Vcs quoted in square brackets hold after 3-family 
unitarity has been imposed. 

The future before a i/MC could contribute can be sketched as follows: 

• Some improvements can be expected over the next several years in direct 
extractions of Vcd and Vcs (i-e. those that do not impose 3-family unitarity 
constraints). However it is very unlikely that they could come close to the 
2 - 3 % level. 
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• A combination of more detailed data on semilcptonic B decays and further 
refinement of heavy quark expansions will yield very significant improve- 
ments in Vcb and V^^; not unreasonable expectations are 

pre—vMC 

~ 10 - 15%. (58) 

Our discussion suggests that potential i/MC analyses could have an essential 
impact on central aspects of the Standard Model by meeting the goals stated 
in the beginning of this section: 

• provide a systematically different determination of \Vcs\ and |V"c6| that is as 
good as can be achieved in charm and beauty decay studies: 

0(few%), lAKbl vMC ~ 

0(few%) ; (59) 

• yield a value for \Vcd\ through direct observation that is about as good as 
otherwise achieved only through imposing 3-family-unitarity 

0{1%) ■ (60) 

• lower the theoretical uncertainty in \Vub\ considerably: 

0(1%) . (61) 

This would enable us to predict various CP asymmetries in B decays with 
order 1% accuracy, thus calibrating the experimental results expected from 
next generation experiments like LHC-B and BTeV and allowing us to ex- 
haust the discovery potential for new physics in B decays; 

• improve dramatically our numerical information on the CKM parameters 
involving top quarks if a very high energy i/MC could be built. 



4 Precision Electroweeik Studies 

4-1 Introduction 

Neutrino scattering is a natural place to study the structure of the weak 
interaction. Historically, it has played an important role in establishing both 
the basic structure of the weak interaction, particularly with the discovery of 
neutral currents [54], and in providing the first precision tests of electroweak 
unification [55]. It is reasonable, as the high intensity neutrino beams at future 
z/MCs offer the promise of a new level of statistics in high energy neutrino 
interactions, that one considers a new generation of experiments to probe the 
weak interaction. 



49 



At the same time, given the ambitious coUider physics programs of LEP I, 
SLD, LEP II and the Tevatron, which will have been completed at the time 
of a z/MC, and given the physics program of the LHC that will be ongoing, 
the goals of such experiments must be correctly focused. 

4.1.1 Knowledge at the Time of uMCs 

Between today and the advent of a i/MC, it is safe to assume that little 
improvement will have been achieved in the amazingly precise measurements 
of decay and production parameters [56]. However, it is very likely that 
LEP II and the Tevatron Run II will have produced a measurement of 
to a 20 MeV/c^ precision[58,59] and that the mass of the top quark will be 
known with a precision of 1 GeV/c^[59]. 

Furthermore, LEP-II has found experimental hints of a possible Higgs boson 
at a mass around 115 GeV/c^ which, although far from being a sure bet, might 
possibly be confirmed at either the Tevatron or the LHC. Alternatively, one 
of these hadron colliders might instead discover the Higgs at a higher mass or 
else something completely different. Of course, the odds of such a discovery 
are much higher if a Higgs exists close to the value hinted at by the LEP-II 
events, as present electroweak fits also seem to suggest [56]. 

An observation of a Standard Model Higgs, along with a precise prediction of 
its mass driven by the W and top quark mass measurements, will make an 
elegant piece de resistance of the electroweak physics program at the energy 
frontier for the period between 1990 and 2010. 

During this time, there may also be improvements in low energy tests of the 
electroweak Standard Model. Atomic parity violation experiments may be able 
to yet again improve significantly in their ability to provide precision tests of 
weak interactions if experiments utilizing trapped unstable Prancium become 
possible[60]. Also, SLAC E-158[61], a proposed polarized MoUer scattering 
experiment, may be able to probe sin^^^K ^ precision of 0.0008 at ~ 
10-2 GeV2. 

4.1.2 Goals of Fixed- Target Electroweak Physics Programs 

The primary goals of a low energy test of a high energy theory remain largely 
the same as they have been in the past. Unification of the interactions of the 
on-shell weak bosons with their low energy manifestations in weak interactions 
at high precision remains an appealing and elegant test of the model. 

The low energy experiments also allow access to some aspects of the theory 
that cannot be readily observed at high energies. For example, a demonstration 
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of the running of the weak couphng strength now only awaits precise low 
energy data [62], given the high precision at the scale of the weak boson masses. 
Another powerful use for low energy data is the sensitivity to interference 
between new physics and tree-level processes. For example, if Z' bosons are 
discovered at the energy frontier then observation of the interference between 
the Standard Model Z and the Z' at low energy may be one of the most 
powerful tests constraining models that relate the two interactions. 



1^.1.3 Electroweak Processes with Neutrinos 

To provide a precision test of the Standard Model, a process must be rea- 
sonably common and precisely calculable. The two useful neutrino interac- 
tion processes for these sorts of studies are measurements of sin^ 9w through 
neutrino- nucleon DIS (neutrino-quark scattering), and neutrino-electron elas- 
tic and quasi-elastic scattering. The former, of course, wins on large cross 
section, but the latter wins on simplicity of target and therefore minimal the- 
oretical uncertainties. The possibilities for using these processes are described 
below. 

Electroweak processes not considered here include neutrino tridents - i.e. three 
lepton final states resulting from internal conversion of a virtual photon - 
which arc effectively a test of boson-boson scattering. While these processes 
are interesting, they do not provide a stable basis on which to form a preci- 
sion probe of the model because of large theoretical uncertainties in the cross 
section. See Sec. 5 for further discussion of this and other processes. 



4-2 Elastic and Quasi- Elastic Neutrino- Electron scattering 



Perhaps the most promising reaction for measuring sin^ 9w at a z/MC is u — e 
scattering. Neutrino-electron scattering possesses one significant advantage 
over DIS for precision electroweak studies, namely that the target is point-like 
and its structure does not introduce uncertainties in extracting the parameters 
of the fundamental interaction from the observed cross sections. 
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4-2.1 Survey of Neutrino- Electron Scattering Processes 



Several u — e scattering reactions will occur in the muon and electron neutrino 
and antineutrino beams at a i/MC: 

v^e~ u^e~, (62) 

u^e" Ue^i~, (63) 

z/^e" u^e~, (64) 

Uee~ i^ee~, (65) 

Pee' PeC", (66) 

PgC" — > z/t-t", ltd . . . (67) 

Reactions 62 and 64 are purely neutral-current processes mediated by the 
exchange of a Z^. Reaction 65 has both charged current {W^ exchange) and 
neutral current components, and reaction 66 has neutral current t-channel 
and charged-current s-channel components. Reactions 63 and 67 can result in 
the production of a single muon and therefore have a significant muon mass 
threshold in the cross section. Figure 8 shows the Feynman diagrams for these 
processes. 

Because of the small ratio of the electron to proton mass, the cross section for 
neutrino-electron scattering is much smaller than that for neutrino-nucleon 
DIS. The leading order differential cross section for neutrino-electron elastic 
scattering with respect to y = Eg/ Ei, is given by 



da , _ _ , G%s 



[ue — > z/e 



dy TT 



9i + gi{l-yY , (68) 



where the center-of-momentum energy, s, is well appoximated by s ~ 2meEi, 
when E^ ^ rric, where terms of OimelEj,) are neglected and where ql and gR 
are process dependent because of their exchange in the neutral current process 
under v ^ v and because of the addition of the charged current process for 
electron neutrino induced reactions. The values of gi and gR are given in 
Table 4. Numerical values for the cross sections after integrating over y are: 



a{ve- ve~) = 1.6 x 10"^' x E^[GeV\ x 



2,^2 

9l + -^Qr 



(69) 



where the values for the final term are given in the final column of Table 4. 
Radiative corrections for this process have been calculated to 1-loop [63], and 
theoretical techniques exist to extend this calculation to higher orders. 

The differential cross section for inverse muon decay, Eq. 63, is 

da Gljs - ml) 

-[y,e - = (70) 
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w 



(C) 





Fig. 8. Feynman diagrams contributing to the — e scattering processes of Eqs. 62 
to 66: (a) NC u-e elastic scattering (for Eqs. 62, 64, 65 and 66) (b) CC i^gC scattering 
(for Eq. 65), (c) CC z/eC annihilation (for equation 66), and (d) inverse muon decay 
(for Eq. 64). 



Reaction 


9L 


9R 


si + 5^1 




-i + sin^ Qw 


sin^ Qy^ 


0.091 




sin^ ^VK 


-\ + sin2 dw 


0.077 


Vf,e~ — > 


\ + sin^ 


sin^ 


0.551 




sin^ 


\ + sin^ 6*1^ 


0.231 



Table 4 

The coupling coefficients, gi, and g'/j, in Eq. 68 for the neutrino-electron scattering 
processes of Eqs. 62, 64, 65 and 66, respectively. The numerical values in the final 
column correspond to the combination of couplings that appears in the tree-level 
total cross section for sin'^ Q^/y = 0.23. 

and the differential cross section for z/gC" — ^ t'^t/U" is 



da 
dy 



Gl{s - ml] 

4:71 



m: 



Ml 



S 



(71) 
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Fig. 9. Integral cross section for neutrino-electron scattering processes above any 
chosen cut on the inelasticity variable, y > ycut, and assuming E^, = 30 GeV. 
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Fig. 10. Change in the integral cross section with respect to sin^ 6w for neu- 
trino-electron scattering processes above y > ycut, assuming E^, = 30 GeV. 

4-2.2 Current Measurements of siri^ 6w from Neutrino-Electron Scattering 

The best measurement of neutrino-electron scattering to date was performed 
in the CHARM II experiment in the CERN SppS neutrino beam. The beam 
was predominantly z/^ and z/^, with an event sample of 2000 events in each of 
the z/^ and z/^ beams. This led to a measurement of the weak mixing angle 
of [64] 



sin^ 9w = 0.2324 ± 0.0058(stat) ± 0.0059(syst). 



(72) 
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Not surprisingly, systematic errors primarily result from normalization and 
background uncertainties. 

4-2.3 Overview of the Measurement Technique at a uMC 

The signature for u — e scattering is a single electron with very low transverse 
momentum with respect to the neutrino beam direction, pt < \fm^El,. There- 
fore, the measured quantity to be converted to a cross section is the number of 
observed events consisting of a forward-going electron track with no hadronic 
activity and with an energy above some defined threshold value, £^cut- 

In order to convert the event count to a cross section, the detector efficiency 
must be determined, backgrounds must be estimated and subtracted, and the 
integrated neutrino flux must measured and/or calculated. 

In the discussion that follows, it will be seen that the physics sensitivity, 
backgrounds and flux normalization procedures will all differ markedly be- 
tween the v^Vi. and v^lV(, beams. Experimental runs with the latter beam will 
have a greater statistical sensitivity to sin^ but the v^Vf. beam will pro- 
vide two experimental advantages: (1) the possibility of flux normalization 
using the muons produced from inverse muon decay and by z/g ~ e~ annihi- 
lation processes; and (2) the background from quasi-elastic electron neutrino 
scattering produces positrons rather than electrons, which can potentially be 
distinguished from the signal electrons by determining their charge sign. 

4.2.4 Statistical Sensitivity 

The number of signal interactions is related to the cross section for the process, 
(j{E^), and to the neutrino flux through the flducial volume of the target, 
^(E^), through 



where the theta function is zero (one) for an argument less than (greater than) 
zero and the proportionality factor is determined by the mass depth of the 
target. 

Figure 9 shows the integral cross section for the reactions above y > ycut, and 
Fig. 10 shows the change in the cross section above y > ycut as a function of 
sin^ 9w (A value of ycut relatively close to zero should hkely be attainable by 
using a dedicated detector for this analysis.) 

The statistical sensitivity to sin^ 9w in any given channel is proportional to 
^/a / {da / ds\T? 9w) 1 and is shown in Fig. 11. In a neutrino beam produced by 




(73) 
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0.5 1.0 

Ycut 

Fig. 11. Statistical uncertainty in sin^ 6w for any chosen value of ycut, from neu- 
trino-electron scattering in beams of either f^, fe, or z^g. An integrated luminosity 
of 10^^ cm~^ at a beam energy of E^, = 30 GeV has been assumed. The values can 
be scaled to other neutrino energies by noting that the measurement's statistical 
uncertainty, for a given integrated neutrino flux through a specified detector, is 
proportional to the inverse square root of the average energy. 

muon decays, the observed rate of visible electrons will include scattering of 
both neutrinos and antineutrinos. 

Note that an undesirable feature of the fi~ beam for measuring sin^ 9w is that 
the dependences of the integral cross sections on sin^ 9]^ for and z/g have 
opposite signs. The resulting sensitivity in integral cross sections for beams 
from /i^ decay is shown in Fig. 12, using the adequate approximation that the 
muon and electron type neutrino fluxes are assumed equal. 

As an aside for a neutrino beam from a polarized muon beam, the exact sensi- 
tivity for this sort of summed measurement would depend on beam polariza- 
tion, particularly for the n~ beam. Due to potential uncertainties in measuring 
the level of any muon beam polarization, better measurements of sin^ 9\y may 
well be obtained from muon beams where the polarization is identically zero 
or can at least be shown to average to zero over the course of a fill. 

4-2.5 Detector Design and Background Rejection 

The best types of detectors for detecting the signal of low pt single electrons 
are likely to be based around kiloton-scale active targets with inherent track- 
ing capabilities and a high-rate capability. In order to minimize the level of 
confusion between the signal process and background events with 7 — > e~e~^, 
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Fig. 12. Same as Fig. 11 except the event statistics has been summed over the 
experimentally indistinguishable contributions from the two neutrino components 
in the "''-induced v^Ue beams and the //"-induced v^De beams produced from muon 
storage rings. 

the target should be composed of only low-Z elements in order to maximize the 
radiation length, and should contain very little dead material. The incorpora- 
tion of a magnetic field to identify the lepton charge would further be helpful 
so as to reduce backgrounds from v^, charged current interactions, and a lep- 
ton charge measurement would also provide a cross-check of sign-symmetric 
detector backgrounds, such as 7 — > e+e~. 



For a detector with all these capabilities, rare low pt backgrounds such as co- 
herent single 7r° production, which were a significant problem in the high-mass 
CHARM II neutrino detector [64], should not be difficult to identify and/or 
subtract on a statistical basis. 



In order to also remove quasi-elastic v^N scattering backgrounds, which can- 
not be separated by electron charge sign identification in the z/^z/g beam, the 
detector will need to be capable of resolving the different pt distributions: z/g 
quasi-elastic scattering off nucleons has a characteristic pt scale of y/rrii^iE^, 
i.e. larger by a factor of ^Jmly/me — 43 than the signal process. In this case, 
signal and background suppression can be achieved by fitting the observed 
single electron pt distribution, which therefore must be measured with a pt 
resolution much better than y/m^Ey. To give a numerical example, the quasi- 
elastic cross section off nucleons at = 30 GeV on an isoscalar target is 
approximately 5 times greater than the inverse muon decay cross section and, 
in this case, a pt cut at 100 MeV would leave a well-measured background of 
about 10% under the inverse muon decay peak. 
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« 100 cm 

Fig. 13. Monte-carlo generated simulation of a high energy neutrino-electron scat- 
tering interaction in liquid helium. The distorted scale is indicated by the 10:1 ratio 
in the relative lengths of the x and z axes. The solid lines are electrons and positrons 
and the dashed lines are photons that would not be seen in the detector. The view 
is perpendicular to both the beam direction and to a 0.1 Tesla magnetic field. See 
text for further details. 



These demanding and specialized requirements suggest a using dedicated de- 
tector rather than, e.g., the general purpose detector of Fig. 2. A natural 
choice is a time projection chamber (TPC) filled with a one or other of several 
candidate low-Z liquids. 



A TPC using the lowest-Z element, liquid hydrogen, may unfortunately be 
ruled out because of insufficient electron mobility, although the possibilities for 
liquid hydrogen TPCs are again attracting some attention [65]. Liquid helium 
also suffers from poor mobility and potentially difficult operation because 
it lacks the ability to self-quench; however, it deserves further consideration 
because it has a radiation length of 7.55 m, allowing very well resolved events. 

Fig. 13 gives an example of a Monte Carlo-generated event in liquid helium. 
The interaction occurs at mid-height at the right-hand side of the figure and, 
in this typical case, the primary electron track travels easily sufficient distance 
to establish its initial vertex, direction, sign and the absence of extra tracks 
emerging from the vertex. Since the event is contained, the primary electron's 
energy could also be cross-checked calimetrically. 

Liquid methane appears to be another good candidate for the TPC medium 
as its favorable electron transport properties have led to it being suggested for 
TPC detectors of up to several kilotons [66]. It is liquid at atmospheric pressure 
between -182.5 and -161.5 degrees centigrade and has a density of 0.717 g/ cm^ 
and a radiation length of 65 cm. Heavier alkanes that are liquid at room 
temperature, such as octane, would be superior for safety and convenience if 
their electron transport properties were found to be acceptable. 

Even if most DIS events are trivial to distinguish from the signal events, 
the detector must still be able to cope with the high interaction rate from 
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the dominant background of DIS neutrino-nucleon interactions, which has 
a cross section three orders of magnitude larger than the signal processes. 
This problem is made even worse for the long drift times typical of the TPC 
geometry because the interactions from up to hundreds of turns (depending 
on experimental details) may pile up in the TPC read-out. This suggests the 
need for additional fast read-out of the events, which could plausibly come 
from, e.g., planes of scintillating fibers within the TPC volume. 

Clearly, any such detector for the neutrino-electron scattering analysis must 
satisfy very stringent experimental requirements and its design and construc- 
tion will be major projects. 

4.2.6 Flux Normalization for Neutrino-Electron Elastic Scattering 

Normalization of the cross section is also a significant issue since this probe 
of weak couplings is only as good as the normalization of the beam fiux. 

Normalization to the muon beam fiux itself is a possibility which would work 
for both z/^z/g and u^Ug beams. Theoretical predictions of the decay process 
would likely not limit this normalization technique. Instead, the ultimate ac- 
curacy should depend on uncertainties in the measured number of muons in 
the ring and on the beam beam dynamics, such as spot sizes, orbits, diver- 
gences and polarization if this doesn't average to zero. It would require a de- 
tailed analysis for a precisely specified muon ring design to determine whether 
or not muon beam measurements and modeling could predict z/MC neutrino 
fiuxes at the 10^"^ level required to be useful for this analysis. 

For the f^Pe beam, an alternative candidate for normalization is single muon 
production in neutrino-electron scattering through the processes of Eqs. 63 
and 67. Like the signal, the absolute cross sections for these normalization 
processes are also extremely accurately predictable but with the crucial dif- 
ference that they don't depend on sin^ 9^. 

It is conceivable that this absolute normalization in the Ufj7,, beam could then 
also be transferred to the u^Ug beam by using the ratio of quasi-clastic events 
to provide a relative normalization between the two beam types. 

4.2.7 Sensitivity to New Physics Processes 

For a beam, the cross section calculations above show that, if the very chal- 
lenging experimental systematic uncertainties can be satisfactorily addressed 
then sensitivities of approximately 6sm^9w ~ 0.0007 would be reached for 
an integrated luminosity of 10^^ cm~^ and a mean neutrino beam energy of 
30 GeV, which corresponds to approximately 1.5 x 10^ DIS charged current 
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events (c.f. Table 2). 



The /x"*" beam's statistical sensitivity would be about a factor of three better, 
with the caveat that beam flux normalization and experimental backgrounds 
are both even more challenging. With an integrated luminosity of 10^^ cm~^, 

expected event sample sizes would be approximately 1.5 x 10^ for a fi~ beam 
and 3 x 10^ for a /j^ beam. Normalization and background uncertainties must 
therefore be kept at the few times 10~^ level in order to achieve this precision. 

Such a measurement could be used to probe for hints of physics beyond the 
Standard Model by, for example, interpreting it in terms of a sensitivity to a 
high mass contact interaction with a Lagrangian of the form 

^= E ^^{^kH.)eH.^^eH^i7H^7f.'^H., (74) 

He{L,R} ^^HeH^ 

where the H indices represent helicity states of the electron and neutrino. 
The statistics given above would probe contact interactions at mass scales 
A ~ 10 TeV, again assuming that statistical uncertainties dominate. 

4-3 sin^ 9w from Deep Inelastic Scattering 

Measurements of sin^ 9w in neutrino-nucleon DIS experiments, using the neu- 
trino beams that have been available from n/K decays, have already provided 
an excellent testing ground for the Standard Model. Nuclcons, however, make 
for a most unappealing and difficult target, and it is necessary to consider 
ratios of observable processes in order to make sense of the results. 

4-3.1 Previous Measurements 

The CCFR u experiment extracted sin^ 9w [67,68] through a measurement of 
the ratio of the cross sections for NC and CC interactions, as expressed in the 
Llewellyn-Smith formula [69] : 

where the value of the parameter p depends on the nature of the Higgs sector 
and has the value p = 1 in the Standard Model. This method, although it 
removed much of the uncertainty due to QCD effects in the target, does leave 
some rather large uncertainties associated with heavy quark production from 
the quark sea of the nucleon target. 
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The successor experiment, NuTeV (FNAL-E815), has improved upon the 
CCFR measurement by using separate neutrino and antincutrino beams. Sep- 
aration of neutrino and antineutrino neutral current events allows the utiliza- 
tion of the Paschos-Wolfenstein relationship [70] : 



'NC ~ ^NC ~ ^-^^ 2 „• 2 , 



'^cc ^cc i r \z 



where 



R~ is, by construction, sensitive only to scattering from valence quarks in the 
proton, and this considerably reduces the theoretical uncertainties associated 
with the target. NuTeV has presented a preliminary result[71] of 

sin^ 9w = 0.2253 ± 0.0019(stat) ± O.OOlO(syst) . (78) 



There were two dominant systematic uncertainties in the CCFR experiment: 
1) Ue flux and 2) CC charm production. These two major systematic uncer- 
tainties in CCFR were reduced in the NtiTcV experiment by utilizing sign- 
selected neutrino beams, leaving event statistics as the dominant remaining 
uncertainty. 

These past neutrino fixed target experiments used dense calorimetric neutrino 
targets in order to increase the interaction rate, and such targets did not 
allow one to distinguish between electron-neutrino-induced charged current 
interactions (CC) and neutral current (NC) interactions. Such experimental 
set-ups would be fatal for i/MC analyses with 2-component i/^Pg and Pf^Ue 
beams and, as will be discussed further; a high performance tracking target 
such as that in Fig. 2 is instead indicated. 



4-3.2 The Experimental Eoctraction of sin"^ 9 w 

As in previous neutrino experiments, the measured quantity used to determine 
sin^ 9w at z/MCs will be a ratio of NC to CC DIS events. However, the NuTeV- 
style ratio of Eq. 76, R~, will not be accessible in the 2-component beams at 
i/MCs because NC events from neutrinos and those from antineutrinos will 
not be distinguishable on an event-by-event basis. Instead, the experimentally 
accessible NC-to-CC event ratios for both the z/^z/e and P^z/g beams essentially 
correspond to linear combinations of the Llewellyn-Smith ratios for neutrinos 
and antineutrinos that were given in equation 75. 
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Table 5 

Estimates of sin^ Q-^ uncertainties in a i/MC analysis, compared to those for the 
NuTeV preliminary result [71] . See the text for details on how the estimated uncer- 
tainties from z/MCs were arrived at. 



The relevant ratios for i/MCs will be: 

a(z/^, NC) + a(l7e, NC) _ + grRp 



Rfx — 



a{v^,CC) + a{v,,CC) 



for the y^^e beam, and 



aiv^, NC) + NC) R, + g-^rR^ 



a{V^,CC) + a{ue,CC) 



1 + g~^r 



(79) 



^0) 



for the u^Ug beam, with r previously defined in Eq. 77 and g the energy- 
weighted flux ratio between Pg and i/fj, in a i/^Pe beam or - equivalently for a 
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non-polarized beam - between and in a v^v^ beam: 



<x>^^ J ^E.)E.dE. J ^E^jEj^dEj,^ 



(81) 



Equations 79 and 80 have made use of lepton universality, which implies that 
UgN scattering cross sections become equal to those for Uj^N at energy scales 
well above the electron and muon masses. The second of the two equations 
differs from the first only in the replacement of g by g~^. 

An analytic calculation [9] gives the value g = Q/7 for the neutrino beam 
produced from an idealized pencil beam of unpolarized muons. It follows that 
the numerical values of the measurements from the i^^Pg and P^i^e beams will 
be nearly identical: 

R^- ~ 0.330; ~ 0.332, (82) 

where we have used the predictions from Eq. 75 of = 0.317 and R'^ = 0.359 
for sin^ 9w = 0.225. Of more experimental relevance, the statistical sensitiv- 
ities to sin^ 9w are also nearly identical, as is indicated by the logarithmic 
derivatives: 

1 dRu- ^ 1 dRn+ , ,„ 

^ - -1.55; , . = -1.47. 



Ra- dsm 9w Ru.+ dsm 9w 



The numerical similarities between the complementary variables i?^- and , 
from i/^Pe and ^nfe beams respectively, mean that the two measurements can 
be regarded as nearly identical from a physics standpoint but with slightly dif- 
ferent experimental systematics due to the approximate interchange of electron 
and muon energy spectra in the CC final states. 

Because of the different kinematics for neutrino vs. antineutrino interactions, 
the CC event sample from the u^V^ beam will contain a softer spectrum of 
primary muons than electrons and vice versa for the z/^z^e beam. The com- 
parison of two theoretically similar measurements with different experimental 
challenges will be a valuable cross-check on the analyses. In this respect, it is 
helpful that muon storage rings are likely [9] to have the capability to reverse 
the polarity of the ring to choose between v^V^ or P^fg beams at any given 
time. 



4.3.3 Detector Requirements 

Any uN DIS measurements of sin^ 9yv at a z/MC would be expected to be a 
large experimental extrapolation from today's measurements. The most de- 
manding requirement on the detector for this analysis will be the ability to 
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efficiently distinguish CC events, with their primary electrons or muons, from 
the purely hadronic events produced in NC interactions. 

The large component of electron (anti)neutrinos in both the V^v^ and v^V^ 
beams at i/MCs rules out use of traditional calorimctric neutrino target / detectors 
since these cannot easily distinguish i/g-induced CC interactions from NC in- 
teractions. 

In contrast, a CCD-based tracking target and general purpose detector such as 
that of Fig. 2 appears to be well suited to the requirements for this analysis be- 
cause of its expected good performance [9] in identifying both primary muons 
and electrons and its further ability to control backgrounds from secondary 
electrons or background muons from pion decays that fake primary leptons. 
Even so, it might even be helpful if the electron identification capabihties of 
such a detector were further bolstered by incorporating transition radiation 
detectors directly downstream from the tracking detectors. 

4.3.4 Estimated Uncertainties 

Table 5 displays a comparison between the uncertainties in the uN DIS mea- 
surements of sin^ Ow from today's most precise measurement [71] and rough 
estimates of the corresponding uncertainties from a uMC measurement, which 
will now be discussed in turn. 

The statistical uncertainty from the roughly one million events at the NuTeV 
experiment was the largest uncertainty in that sin^ 9w analysis. Table 2 sug- 
gests vMC event statistics of 10^ events or more, corresponding to a reduction 
in statistical uncertainty by at least a factor of 30 and pushing the absolute 
statistical uncertainty in sin^ 9w^o below 0.0001. 

Turning now to experimental uncertainties, the NuTeV uncertainty of 0.0004 
due to Ve flux was relevant only for calorimctric neutrino targets and will no 
longer exist for the tracking targct/dctcctors discussed for z/MCs that will 
be capable of distinguishing, on an event- by-event basis, between NC inter- 
actions and the CC interactions of fg's. Uncertainties from energy scale and 
event length in the NuTeV analysis will also be irrelevant for vMC detectors 
because they were associated with NuTeV's statistical event-length method of 
separating NC from CC events. 

The improved vMC method of identifying primary leptons on an event-by- 
event basis will instead have to contend with uncertainties in the identification 

efficiencies for the primary leptons that distinguish CC from NC events. Every 
misidentification moves that event between the numerator and denominator 
of the experimental ratios of Eqs. 79 and 80, so it is clear that the fractional 
uncertainty in the level of misidentifications must be reduced to well below 
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the 10""^ magnitude desired for the fractional uncertainty in sin^^^y. Some 
confidence that this might be achievable comes from the very high lepton 
identification efficiencies for the detector scenario discussed above. 

Both the rejection of backgrounds and the positive identification of the pri- 
mary lepton are generally more difficult for low energy leptons, so the sin^ 9w 
measurement would benefit from using cuts on the minimum lepton energy. 
The value of this energy cut must be balanced against increasing theoreti- 
cal uncertainties as progressively more of the event sample is cut away. In 
any case, estimation of the identification efficiency for primary leptons may 
well be the largest experimental uncertainty in a measurement dominated by 
theoretical uncertainties. 

An improved understanding of several potential theoretical uncertainties will 
be required to attain a sin^ 9iy measurement that could be meaningfully in- 
terpreted as equivalent to a sub— 25 MeV W mass uncertainty. In particular, 
calculations and/or measurements to minimize the charm production uncer- 
tainty, higher twist effects, radiative corrections, and uncertainties in the lon- 
gitudinal structure function will need to be dealt with. 

In this respect, it is helpful that several of the theoretical uncertainties can 
be calibrated using the same z/MC event sample as used for the sin^ 9w anal- 
ysis. Hence, the enormous increase in statistical power of z^MCs over today's 
neutrino experiments should also help to minimize some of the systematic 
uncertainties in Table 5. 

Good examples of theoretical uncertainties that are amenable to experimental 
calibration arc the large uncertainty due to charm mass effects and the related 
uncertainties in estimating the charm and strange seas of the nucleons. A Next- 
to- Leading- Order (NLO) fit to the charm mass from CCFR [17] obtained rric — 
1.71 ± 0.19 ± 0.02 GeV/c^, corresponding to a charm production uncertainty 
in the CCFR sin^ 9w measurement of 6 sin^ 9w = 0.003. Since the statistics at 
a z/MC might be three or more orders of magnitude larger than in the CCFR 
experiment. Table 5 somewhat arbitrarily chooses an improvement by a factor 
of 10 on the CCFR uncertainty. This is less improvement than the factor of 
30 or more that would be predicted from straightforward statistical scaling 
but a careful analysis would be required to establish the actual level at which 
residual theoretical uncertainties set in that cannot be calibrated away using 
the experimental data. 

Strange quark mass effects in i/iu —>■ i~s and u^s — > provide a much smaller 
theoretical effect that fell below the uncertainty threshold for the NuTeV 
analysis but whose corrections and uncertainties would need to be checked 
for a z/MC analysis. Although presumably a very small correction, its effects 
are difficult to reliably establish because lattice gauge calculations predict that 
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nis has an awkward value: rus ^ 300 MeV~ h-qcD- 

Higher twist effects are assessed as one of the larger theoretical uncertainties 
in today's measurements. However, a recent study [72] indicates that most of 
the higher twist effects might be able to be reinterpreted as higher order QCD 
corrections that can be determined from the structure functions measured 
in the z/MC data sample. For the theoretical precision required at a z/MC 
measurement, it may be necessary to evaluate and correct for the residual 
small effects from "radiative higher twist processes" such as v^n — > '^^~p. The 
radiative photon present in CC events generates a CC/NC asymmetry and can 
boost the apparent of events to high enough values to evade cuts designed 
to suppress higher twist corrections. 

As was already mentioned in Sec. 2, cancellations of theoretical uncertainties 
by applying isospin invariance relations are very important for reducing the 

uncertainties in sin^ 6y/ and in other analyses in uN DIS experiments and this 
is the motivation for using neutrino targets that are approximately isoscalar. 
This theoretical handle was useful for both the CCFR and NuTeV analyses, 
using a neutrino target detector with a neutron excess of {N — Z) / {N + Z) ~ 
0.0567 ± 0.0005. Isospin invariance relations should be even more applicable 
for, e.g., a CCD target at a z/MC since the silicon substrate of CCD's has a 
neutron excess of only 0.3% and even the residual target components will have 
small non-isoscalarities: aluminum conductor has a neutron excess of 3.6% and 
the dominant carbon component of the support structure has an excess of less 
than 0.1%. 

Indeed, the z/MC neutrino target might be sufficiently isoscalar that uncertain- 
ties due to the neutron excess might fall below those due to nucleon isospin- 
violating effects [73] arising from electromagnetic and quark- mass contribu- 
tions that break the generally assumed isospin invariance relations between 
the up and down quarks in protons and neutrons: Up{x, Q^) = dn {x, Q^) and 
Un{x,Q'^) = dp {x,Q^). Such effects will be present even in a deuterium tar- 
get [9] or other pure isoscalar targets. 

The theoretical uncertainty due to electroweak radiative corrections in the 
NuTeV experiment, Ssin^Ow — 0.00051, deserves further theoretical atten- 
tion for a z/MC since it does not appear to be amenable to data-based reduc- 
tions from the improved statistics and experimental conditions at a z/MC and 
so it might well dominate the total uncertainty at a z/MC if the theoretical 
calculations are not improved upon, as will now be discussed. 

The current state of electroweak radiative corrections suffers from the fact that 
only one attempt at a complete calculation exists, by Bardin and Dokuchaeva[74], 
which, curiously, is unpublished. This calculation includes electroweak effects 
to one loop, is lengthy and complicated, and it would benefit from independent 
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confirmation. Further, it has several theoretical shortcomings, as follows. The 
Bardin-Dokuchaeva result depends explicitly on quark masses, which intro- 
duces a spurious model dependence. Similar calculations of EW radiative cor- 
rections in W mass production in pp annihilation show that quark mass effects 
can be absorbed into parton distributions and fragmentation functions [75]. 
The calculation is only approximately valid to leading log in QCD. In particu- 
lar, it neglects scaling violation effects in important diagrams involving muon 
bremsstrahlung. It also neglects contributions from longitudinal partons and 
effects from the target mass. Further, it does not incorporate heavy quark 
effects in CC and NC charm production and, finally, it is not differential in 
the final state radiated photon momentum vector. 

A i/MC measurement of sin^ 9w would benefit greatly from a new EW radiative 
correction program with the aim of reducing the residual error on the effective 
W mass to ±1 MeV. To reach this level, it may be necessary to re-sum large 
lepton logs, to include order aas contributions, and to apply EM radiative 
corrections to input parton distribution function sets. Such a program should 
also provide the capability for the explicit generation of jiiq final states. 

4-3.5 Comparisons of Expected Precisions from sin^ 9^ Measurements in Dif- 
ferent Experimental Processes 

To summarize the content of Table 5, the uncertainty on a i/MC DIS mea- 
surement of sin^ 9w might well be equivalent to on the order of a 25 MeV 
uncertainty on the W mass and could perhaps improve on this if the theoret- 
ical uncertainties due to radiative corrections can be controlled. 

In order to put this in the context of collider measurements at the energy fron- 
tier. Fig. 14 adds the inferred W mass information from a vMC determination 
of sin^ 9w to a plot showing the expected status of coUider W and top quark 

mass measurements by the year 2010. Although the level of the bands and the 
actual slope of the curves might change depending on other parameters in the 
measurements, this figure gives a good idea of the level of accuracy one might 
expect. 

The most precise measurement of My[/ by the year 2010 is expected to come 
from direct measurements at TeV33. The contour represents the 68% confi- 
dence level from the TeV33 expectations of 6M]y = 20 ~ 30 MeV/c^ and 
5Mt = 1 ~ 2 GeV/c^, with J Cdt = 10 fb"^ and using the traditional Mt 
method [76]. As can be seen in Fig. 14, since the errors from both the direct 
measurements and the z/MC are going to be extremely small, and the z/MC 
measurement provides the Standard Model-based band in Myy — plane, 
the measurements would be complementary to each other in testing the Stan- 
dard Model by providing an independent prediction of the Standard Model 
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vs in 2010 




Fig. 14. Expected future restrictions in the size of the allowed regions in a plot of 
W mass vs. top quark mass, following experiments at TeV 33 and at a z/MC, as 
illustrated by the thicknesses of the bands and ellipse. The exact positions of the 
shaded regions are for illustrative purposes only. 

Higgs mass to better than ~ 20%. 

4-4 Summary on sin^ ^vi^ Measurements at uMCs 

This section has demonstrated that precision tests of the electroweak section 
of the Standard Model can be expected to play an important part in the 
physics program of z/MCs. Two possible tests were discussed, one each from 
neutrino-nucleon scattering and neutrino-electron scattering. 

The former process will allow tests of unparalleled statistical accuracy but 
will likely suffer from substantial QCD uncertainties and, perhaps, also from 
experimental systematic uncertainties. By contrast, neutrino-electron scatter- 
ing may be a statistical challenge even at a z/MC, but it offers the possibility 
of a very clean measurement if experimental systematics associated with nor- 
malization and backgrounds can be well enough controlled by an appropriate 
design of the detector and analysis. 
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5 Hare and Exotic Processes 



Despite the impressive direct and indirect searches for new physics available at 
the higher center of mass energies at the Tevatron, LEP, HERA and LHC col- 
liders, some searches remain unique to neutrino experiments and these could 
be improved quite significantly with the much higher event statistics and im- 
proved experimental conditions promised at z/MCs. This section presents sev- 
eral examples of such processes, emphasizing the complementarity of these 
studies to already existing programs. 

Additionally, two rare processes, u^e' annihilation and W/Z-j scattering, 
present novel tests of low energy features of the Standard Model. 

5.1 New Physics Sensitivity 

5.1.1 Flavor Changing Neutral Currents 

Besides testing current predictions of the Standard Model, a ulsAC offers op- 
portunities to search for new phenomena in yet unexplored physical regions. 
An example is neutral current production of a single heavy quark: 

v^N^u^hXB=o. u^N^vJ)Xb=o. (83) 

which would signal the presence of flavor changing neutral current (FCNC) 
processes. These reactions can provide important constraints on new physics 
as they occur in the Standard Model only at the one loop level while new 
physics effects can occur at both the tree level and one loop level. Examples 
of the former include new intermediate bosons with FCNC quark couplings 
while examples of the latter include a wide class of new physics models such as 
supersymmetry and technicolor. Some of these models are already constrained 
from other measurements [84]. Unfortunately, gains in sensitivity to this type 
of new physics increase only slowly with event statistics. For instance, in the 
models with new tree-level FCNC interactions, such as string-inspired models 
with neutral Z' bosons, the FCNC rate due to Z' exchange is proportional to 
l/M^i and so sensitivity to Mz' only improves at best as the fourth root of 
rate increases. 

The FCNC vertices vvsb that contribute to equation (83) will be extensively 
studied in exclusive and inclusive B decays at the B factories as well as at 
CERN [83]. ALEPH has obtained the best bound so far: BR{B vuX) < 
7.7 X 10~^. Within the Standard Model one predicts a value of ~ 4 x 10~^. 
In contrast, studies of the uudb vertex will be extremely challenging even at 
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future B factories. Processes that involve neutrinos, such as 5 — > X^i/P and 
i/^N v^hXB=Qi have the considerable advantage over the corresponding 
FCNC processes involving charged leptons that their rates are mostly deter- 
mined by short-distance physics, which ensures the robustness of the theoret- 
ical predictions. 

, in the event of any observed signal, would 
for the irreducible Standard Model background. 

Single h quark production is enhanced in the Standard Model through the 
GIM mechanism acting with high mass intermediate top quark states: 

a (uN ^ ubXB=o) oc GUmt/M^) \n{mt/M^) (m^jM^) [\Vuf D + S] r] (ml/2ME) , 

(84) 

(with next-to-leading order QCD corrections available [84]) where D and S 
represent the down and strange quark contributions from the nucleon, respec- 
tively, r]{'ml/2ME) represents a kinematic threshold suppression from the 
heavy b quark mass, and the contributions from D and 5* arc likely to be simi- 
lar. Even so, Eq. 84 predicts the FCNC with a uudb vertex to occur only at the 
level of 10~^ of the event sample even for uMCs at high enough energies for 
the B threshold effects to become small. Most likely, therefore, the Standard 
Model backgrounds will instead come from other processes that have been 
experimentally misidentified and goal of the analyses will be to search for new 
physics effects that enhance the FCNC event sample to considerably above 
the predicted background level. 

For single charm quark production, the irreducible Standard Model back- 
grounds will almost certainly be negligible since, besides the CKM factors, 
the production amplitude is suppressed in the Standard Model by m^m^ rel- 
ative to single b production, although with the caveat that the prediction is 
more sensitive to long-distance QCD effects that are not currently calculable 
from first principles. Again, experimentally misidentified events will dominate 
the backgrounds. 

A high performance detector with excellent vertex tagging, lepton identifica- 
tion and reconstruction of event kinematics will be required to cope with large 
background levels involving both CC and NC allowed production of charm or 
beauty. The allowed CC channels, 

Ui,N fjL'cX, (85) 
PeN e+bX, (86) 

will be most dangerous as |/ — > 1 and so the very low energy final state muon or 
electron can escape detection. This background could be suppressed to some 
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extent by imposing a cut on the minimum allowed transverse momentum, pt, 
in the event, which can be large for the signal NC events when the neutrino 
has a large pt but should be zero within the detector resolution for charged 
current events. Neutral current production of heavy quark pairs 



forms the other background when one or other of the heavy quark decays is 
not picked up by the vertex detector because it occurs too quickly or goes into 
an unfavorable final state. This background in particular makes setting any 
stringent limits from FCNC production of charm difficult at z/MCs, even with 
an excellent vertexing geometry such as that shown in Fig. 3. 

The situation is more promising with B decays, since almost every B decay 
gives two chances for detection: at the primary B decay vertex and at the 
decay of the daughter charm hadron. Also, charged current b production may 
be accompanied by a c, and this information can be used as well. 

Even for FCNC B production at a z/MC, the only interesting limit may well 
be for the subset of FCNC B production that occurs off a valence d quark, 
i.e. involving the FCNC vertices uudb. This restriction provides two important 
additional experimental handles: (1) kinematically, almost all of the NC back- 
ground events will be at relatively low Bjorken x while the PDF for d valence 
quarks extends to high Bjorken x, and (2) d valence quarks will produce B~ 
mesons approximately half the time but never a B~^, while the S's forming 
the NC background are sign symmetric. 

Further discussion on the experimental and theoretical issues for uudb FCNC 
searches at uMCs is given elsewhere[14]. 



5.1.2 Generic Four-Fermion Operators 

Ncutrino-nucleon processes at low momentum transfer are sensitive to generic 
four-fcrmion contact terms produced by the high energy neutral current inter- 
actions. Other, fiavor changing, couplings are well constrained by the limits 
on processes like tt — > eu. These four-fermion interactions can be generated in 
a variety of new physics scenarios. Examples include, but are not limited to, 
super symmetric theories with 7^-parity non-conservation, new vector bosons 
Z' which appear in many superstring-motivated models, models with TeV- 
scale gravity, and quark compositeness models. For instance, the low energy 
remnant of a generic high energy electron-quark neutral current interaction 



^t^N U/^bbX, 



(87) 
(88) 
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can be represented by non-zero coupling constants, rj, in the Lagrangian: 

^Nc = \^LL {eLli^ei) {qEYQL) + Vrr {eR'-fi^en) {WI^Qr) 

Wlr {^Iu-^l) {iRl^qR) + Vrl (eRlneR) {qLl^QL)] ■ (89) 

A similar equation can be written for direct neutrino-quark interactions. One 
can use SU (2) symmetry to relate v and e couplings: 

'ILL — 'ILL 1 
'ILR — 'ILR ) 

Vli = vI'r , (90) 

so that i/N interactions can be used to constrain the ry's in the Lagrangian of 
Eq. 89. 

A particular example of a high-energy model that leads to a low-energy La- 
grangian of this type is provided by R-parity- violating ( ^) SUSY, with the 
Lagrangian: 

-f>i44 - did^ui - 4*4'4) + h.c. (91) 

At low values of transferred momenta, this Lagrangian can be approximated 
in terms of local four-fermion interactions: 

J^ed = — 2 — ^LdRdReL H ^ — ^hdRdRVL 

mzi m~i 



Assuming that the squarks of first two generations are degenerate and impos- 
ing SU (2) symmetry constraints gives 



ed _ (Kji? _ (Kji? _ ud /qq^ 



4 X 



Indeed, the best constraint on this coupling, t]!*^ < 0.07lo;|4 already comes 
from the analysis of neutrino-nucleon scattering experiments [85] . Data from 
vMCs should complement new constraints on the new physics contributions 
that will become available from new Tevatron, LHC and DESY analyses. 
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5.1.3 Heavy Neutral Lepton Mixing 



Another opportunity for vMCs to significantly improve on already existing 
bounds on new physics from neutrino experiments is provided by the ability to 
search for the existence of neutral heavy leptons [86]. In several models [87,88], 
neutral heavy leptons are considered heavy isosinglets that interact and decay 
by mixing with their lighter neutrino counterparts. The high intensity neutrino 
beams created by i/MCs provide an ideal setting to search for neutral heavy 
leptons with a mass below 100 MeV. 

It is postulated that neutral heavy leptons, Lq, could be produced from muon 
decays when one of the neutrinos from the decay mixes with its heavy, iso- 
singlet partner. The expression for the number of neutral heavy leptons pro- 
duced in a muon beam is: 

Nl, = TV, X Bri^^-^^) X e x e"^/^- 
fi we 

X ^,( -^0 - detectable ^ ^ _ ^^^^ 
LiQ — > total 

Here is the number of neutrinos produced from muon decays, i?r(/x — > 
Love/ ^ — > /iz/e) is the branching ratio of muons decaying into neutral heavy 
leptons versus ordinary muon decays, L is the distance from the beam-line 
to the detector, 51 is the length of the detector, e is the combined detec- 
tor efficiency and geometric efficiency, r is the Lq lifetime, and Br{Lo — > 
detectable/Lo total) is the branching ratio for the neutral heavy lepton 
decaying via a detectable channel (mainly, Lq vee). 

Note here that the muon has two possible ways of producing Lq's: 

At" ^ -Lo + I^e + eT (95) 
jjT ^ v^^ Lq^ e~ . (96) 

The branching ratio for each of these reactions is given by: 

Br{ii LoV^e/ii I'^ii'ee) = \U^Lf{l - 8a;^ + 8a;^ - a;^ + 12a;^lna;^) 

(97) 

Br{ii f^LQe/n ly^i^ee) = \UeLf{l- Sx^ + Sx^ - + 12a;^lna;^). 

(98) 

Here Xm is defined as mL^^/m^ and |f/(^,e)Lp is defined as the coupling constant 
between the specific type of neutrino and the neutral heavy lepton. 

Once produced, a neutral heavy lepton of such low mass will either decay via 
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three neutrinos (Lq — > lyi/u), or through the channels: 



(99) 
(100) 



The first decay mode involves only the neutral current, whereas the second con- 
tains a mixture of both neutral and charged current interactions. The branch- 
ing ratios for decay processes (99) and (100) have been calculated [8 9]. Note 
that the number of Lq's detectable at the i/MC depends roughly on U^. 

Using the above model with some additional assumptions, one can estimate 
the number of neutral heavy leptons produced at the vMC that then decay 
to two electrons and a neutrino and are detected in an experiment. The plots 
in Fig. 15 show limits for the coupling constants at a z/MC as a function 
of the Lq mass and for a number of different muon energies and detector 
distances. All the plots assume a pure, unpolarized muon beam containing 
10^° muons/year with straight sections such that 25 percent of the muons will 
decay to neutrinos pointing towards the detector. The detector parameters are 
based on the detector for the NuTeV Lq searches [90-92]: 3 meters in diameter 
and 30 meters in length and with enough resolution to detect the e+e~ vertex 
from the Lq decay. Finally, it is assumed that backgrounds are negligible. 

It is seen from Fig. 15 that one achieves the best limits from using relatively low 
energy muon beams. This is a significant improvement over previous neutral 
heavy lepton searches, where limits fail to reach below 6.0 x 10~^ in the low 
mass region. 

The i/MC may prove to be an ideal location to continue the search for neutral 
heavy leptons. The high intensity neutrino beam allows for a neutral heavy 
lepton search to be sensitive to coupling constants in the low mass region. 
In addition, such a neutral heavy lepton program could easily interface with 
an already existing detector utilizing the neutrino beam. It is also clear that 
a neutral heavy lepton search would receive the most benefit at lower muon 
energies, and thus would yield best results at lower energy i/MCs. 



5.2 Studies of Low Energy QCD 
5.2.1 Dee~ Annihilation 

Section 4.2 showed that neutrino-electron elastic and quasi-elastic scattering 
can be used to extract the electroweak parameter sin^ Ow from purely leptonic 
interactions. This section instead discusses the neutrino-electron annihilation 
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Fig. 15. Plots for limits on the (from top to bottom) coupling values \UeL\^, IU/^lUclI 
and {UulI"^ as a function of Lq mass for one year of running under the assumptions 
given in the text. The solid, dashed, and dotted lines represent uMC energies of 10, 
20, and 250 GcV respectively. The plots on the left show limits for L = 12 m, and 
the plots on the right show limits for L = 1 km. 

processes of Eq. 67, 

i>ee~ ^ v^H^ ,UrT~ ,ud . . . (101) 



Electron antineutrino-clcctron annihilation is expected to show considerable 
complexity once the center of mass energy exceeds the threshold for annihila- 
tion into hadronic final states. The process can be compared to e+e~ — >hadrons 
but with the p, and other low- lying vector resonances replaced by the 
n~, p~ , and aj~. The axial component of the weak current produces coupling to 
axial vector resonances and, at low energy, to the pion. These latter couplings 
are absent in electron-positron annihilation and thus the weak annihilation 
of the electron offers a novel complement to electron-positron physics. Direct 
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Table 6 

Threshold neutrino energies for various hadronic final states. 

measurements of Vf.e" annihilation would complement the detailed studies of 
hadronic tau decays that currently provide some of the most powerful tests of 
QCD at low momentum transfer. 

These measurements bridge two asymptotic limits of QCD: the perturbative 
regime operative at high and the chiral limit (to^ = irid = mg = 0) that is 
approached at low ^/s, and would provide some important information about 
the QCD spectral functions used in reducing the theoretical hadronic uncer- 
tainties in studies of CP -violation in kaon decays, particularly in the interpre- 
tation of e'/e measurements [77,78]. Unfortunately, since the t'i^t threshold 
occurs at E,^ — 3.1 TeV, there is a long way to go in neutrino energy in order 
to completely overlap the ^/s region probed by tau decays. 

Table 6 summarizes neutrino energy thresholds for s-channel final states. The 
thresholds for all channels are quite high. Nevertheless, one can see that a 
50 GeV muon storage ring could provide access to the lowest lying final states 
and explore the interesting region near ^/s = rrij^ to provide a clean test of 
PCAC, and a 250 GeV muon beam would extend the reach up to the threshold 
for kaon production. 



5.2.2 W*/Z* -photon Scattering 

Just as intense neutrino beams open up the possibility of directly exploring 
weak annihilation of leptons in analogy to e+e~ annihilation, so too does one 
acquire access to the analog of two-photon physics: the scattering of virtual 
W or Z beams from quasi-real photons in the Coulomb field of the nucleus. 
The general reaction is of the form 

ueA ^ feFfA, 
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where F, / arc fcrmions, and fe = ue or £~ . For Ff = z/^^"*" both CC and NC 
diagrams contribute. Otherwise fi — i~ proceeds through Wj scattering and 
= vi through Z7 scattering. 

The creation of lepton pairs in the Coulomb field of the nucleus is sensitive to 
W — Z interference and provides a direct test of the Standard Model. Previ- 
ous experimental observations are ambiguous[79,80,53]. The purely coherent 
part of the cross section for vi^A i\Vi^i\ at asymptotically large energies 
is [81,82] 

a,,,, (E^oc)^ ^^, EK,,,, (E, A) (1 - x) , 

Ke,e, {E, A) = log ' + - log ' - R,,,, . (102) 

Ke^£2 {^1 ^) is a reaction and nucleus-dependent form factor with j3 {A) ~ 

fj-Q-i — gA"^/^ GeV set by the nuclear size, and p^j^2)P^i^2) -^^1^2 
simple functions of either /3 [A) or final state fermion masses. The factor Z'^a^ 
reflects the coherent electromagnetic nature of the process; and the scale of 
the cross section is set by G'j^PiE. Since the effective center-of-mass energy 
is ^/s ~ y/2jJjE and f3i ^ mg, all possible leptonic final states are accessible 
(including those with r^) to a neutrino beam derived from a 250 GeV muon 
beam. 

The factor x incorporates effects of neutral currents, including interference 
in the £'^£~ final states. A nice electroweak test is to measure x through the 
ratios 

a^^ {v^A fi-fi+u^A) ^ Kf,f, {E, A) 

a,, {u^A ^ ii-e+v,A) ^ ^> K^, {E, A) ' 

where K^f^ [E, A) / K^^e [E, A) will depend only weakly on nuclear form factor 
and energy. 

Hadronic resonances are also possible. We expect similar states to those pro- 
duced in 77 collisions, namely 0""*", 0"*""*", 2"*""*", ...with 7 = and I — 1. Some 
Cabibbo-allowed examples include 

UeA e'alA (103) 
i/eA e~D'^A. 

Unfortunately, some of these hadronic resonances do not provide a unique ex- 
perimental signature. Single pion production, for example, can occur through 
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the diffractive process UgA e n'^A mediated by a collision of the pion com- 
ponent of the virtual W boson with a pomeron from the nucleus or nucleon. 

Experimental backgrounds that have to be controlled for these studies include 
coherent meson production mediated by VF/Z-pomeron scattering. Especially 
tricky are the diffractive D*/Ds channels where the Dg undergoes two-body 

r decay, 

ueA rDt*/D+A r (7) r+ur ^ (7) i'^i^ri^ri^e'- (104) 

Also of concern are backgrounds from inclusive charged current charm pro- 
duction, i/^N —>■ iJ,~cX, where the charm quark fragments into a charmed 
hadron which takes nearly all of the event's hadronic energy and then decays 
semi-leptonically into a final state where leptons carry nearly all the energy. 

5.3 Conclusions on Rare and Exotic Processes at uMCs 

Rare processes that could be studied at a i/MC would probe, both directly 
and indirectly, an energy range from fractions of a GeV to above the TeV 
scale. This would provide important information complementary to the exist- 
ing results in some areas (e.g., in low energy QCD studies, FCNC, contact 
interactions) and could substantially improve current bounds on the parame- 
ters of some new physics models (e.g. heavy neutral lepton searches). 



6 Chcirm Decay Physics 

6. 1 Introduction 

A i/MC will constitute a rather impressive charm factory. Figures 4 and 5 of 
Sec. 1 show that one can expect between 2x10^ and 2x10^ well-reconstructed 
charm events in a total event sample of 10^° events, depending on the vMC 
energy. Several species of charmed hadrons should be produced, with mea- 
sured [93] relative production fractions for the more common charmed hadrons 
of: 

D° : : : AJ = 0.60 : 0.20 : 0.10 : 0.10. (105) 

Also, the and are expected[94] to have comparable production cross 
sections to Aj. The charmed-strange baryons and should be produced 
at levels down by a factor of a few and Jlc should be still less common. 
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The ratios of Eq. 105 arc relatively independent of the neutrino energy for en- 
ergies above 10 GeV and are for production from neutrinos; the correspond- 
ing antiparticles containing anti-charm will be produced from antineutrinos 
in similar ratios, although with differences in the absolute cross sections and 
kinematic distributions. The large asymmetry between and D'^ production 
is due to the prevalence of D* production with its preference for decays into 

As well as providing good all-around event reconstruction, i/MCs will have two 
other distinct and important experimental advantages over all other types of 
charm facilities. Firstly, reconstruction of the charm decay vertex should be 
superior to that at any collider experiment, particularly for the reconstruction 
of the challenging 1-prong charm decays, as was illustrated by Fig. 3. Secondly, 
a uniquely pure and efficient tag of whether the production flavor is charm 
or anti-charm is provided by the 100% correlated sign of the primary lepton 
from the interaction: 

PiQ e+c. (106) 

This section discusses several areas for charm decay physics at uMCs where 
these experimental capabilities should be important. The theoretical interest 
of each measurement will be discussed, and brief summaries of the expected 
experimental techniques and sensitivities at uMCs will be included. Expected 
relative strengths and weaknesses of zvMCs compared to other future charm 
facilities will also be touched on. However, detailed numerical predictions for 
measurement precisions await more extensive feasibility studies than have been 
performed for this report. 

There are also possibilities for B decay physics using the neutrinos from multi- 
TeV muon colliders[14,9], where bb production in neutral current interactions 
should be at the level of 10~^ of the total cross section. Associated production 

of be and cb can also be studied, however the relevant production cross-sections 
are suppressed by approximately two orders of magnitude compared to the bb 
production cross section. 

6.2 Theoretical Motivation for Charm Physics 

It is clear that, from the point of view of Standard Model electroweak physics, 

charm decays represent a decidedly dull affair. First, the relevant CKM param- 
eters are reasonably well known, for the smallness of \Vcb\ and \ Vub\ constrains 
Vcs and Vcd very tightly through three-family unitarity (see Sec. 3). Second, 
]jo _ £)0 oscillations proceed slowly. Third, CP asymmetries are small due to 



79 



the fact that both decaying and final state particles contain quarks of only 
the first two generations. Finally, rare charm decay rates are tiny and, again, 
are dominated by long-distance effects. 

These apparent vices can, however, be turned into virtues. Since the weak 
dynamics apparently hold no secrets, one can employ charm decays as a lab- 
oratory to study QCD in the interface of perturbative and non-perturbative 
dynamics. Also, precisely because the Standard Model promises us no drama in 
charm decays, one can conduct searches for — D^ oscillations, CP violation 
and rare charm decays as probes for new physics with almost no background 
from the Standard Model. 

6.3 Probing Strong Interactions through Charm Decays 

Improved measurements of charm decays are needed for phenomenological 
and theoretical reasons even in the absence of new physics, for the following 
reasons: 

• to improve the data base needed for analyzing B decays one needs more 
precise measurements of the absolute branching ratios of charm hadrons; 

• measurements of the leptonic decay rates Di^g) are required for deter- 
mining the meson decay constants; these decay constants give us quantita- 
tive insight into the dynamics of heavy-light bound state systems and can 
be used for tuning the lattice QCD methods and a more reliable evaluation 
of 5° - S° oscillations; 

• more precise studies of inclusive semileptonic D, Dg, Ac, etc. decays would 
provide us with valuable novel insights into the inner workings of QCD and 
at the same time sharpen our tools for a quantitative treatment of B decays. 

6.3.1 Absolute Charm Branching Ratios 

As the discussion about the charm content in the final state of B decays illus- 
trates, a significant bottleneck in the detailed analysis of beauty decays of the 
b ^ c type is currently caused by the uncertainties in the absolute branching 
ratios of charm hadron decays to specified final states, in particular of Dg, Ac 
and Sc. A z/MC should be well suited to obtaining these branching ratios, as 
we now discuss. This information will be useful even if it is obtained only after 
the next generation of B experiments have accumulated their samples. 

As what is typically the less difficult part of the measurements, the expected 
excellent particle identification and event reconstruction at z/MCs should give 
good capabilities for determining relative branching ratios for each hadron. 
The more difficult task of obtaining the production normalization factors to 
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convert these to absolute branching ratios should then be achieved by fitting 
the experimental decay length distributions in a procedure that was studied 
for the COSMOS (E803) neutrino experiment at Fermilab. 

The COSMOS technique[95] envisions fitting normalization factors to the sev- 
eral known decay exponentials - one for each charmed hadron species - in the 
observed neutral and charged distributions for the variable x = d/p, with d 
the charmed hadron distance to the decay vertex and p its reconstructed mo- 
mentum. It is helpful that the exponential decay constants in this variable are 
well separated for both the charged and neutral hadron distributions: 

170 /im/(GeV/c) 
71 /im/(GeV/c) 
43 /im/(GeV/c) 

27/xm/(GeV/c), (107) 



67/xm/(GeV/c) 
12 ^m/(GeV/c) 

7/xm/(GeV/c) (108) 



and 
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x(S+) 
a;(A+) 



x(D' 



x(^ 



c) 



x{nl 



c) 



Auxiliary information for the fit will be available from particle identification 
in the detector. In particular, the presence of a proton in the final state will 
reliably indicate the decay of a baryon rather than a meson. 

To test the method for the COSMOS environment, exponential fits were per- 
formed [95] for simulated decay length distributions from approximately 14 000 
reconstructed D"*", D'g and charm decays. The fitted statistical uncertain- 
ties for the three species were 3.4%, 12% and 5.4%, respectively. These sim- 
ulations show that statistical uncertainties would be negligible for such a fit 
at a i/MC, which would have several orders of magnitude more events. The 
uncertainties in the charm production rates would instead be dominated by 
uncertainties in modeling the level of vertexing inefficiencies. Hopefully, these 
uncertainties could also be made small due to the favorable vertexing geome- 
try shown in Fig. 3 and to the considerable potential for using the data itself 
to estimate the inefficiencies. 

Another area where z/MCs can be expected to make significant or even unique 
contributions is in the analysis of final states that contain more than one 
neutral hadron, e.g., 

D° ^ tt+tt-ttV ; D+ tt+ttV ; ^ 7r+7r°77 etc., . (109) 
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Even all neutral final states like 

^ 27r°, 37r° (110) 

might become observable. 

Such neutral-rich channels are rather elusive for the usual e^e~ annihilation 
and photoproduction experiments. A i/MC could access these modes through 
the expected sample of 10^~^ NC-produced cc events. Vertex tagging one of 
the charmed hadrons would allow a search for such decay modes in the other. 
Filling in these 'white spots' in the map of charm decays would close or at least 
narrow the gap between exclusive and inclusive decays and thus can provide us 
with important lessons on how quark-hadron duality is realized in subclasses 
of total decays. For example: a quark based description leads to the prediction 
that the (Cabibbo suppressed) inclusive rates driven by c ^ ssu and c ddu 
should practically coincide since m^, uis <^ nic. Yet exclusive channels like 
K^K'' and — > 7r"'"7r~ do not at all follow this expectation! Duality 
suggests that a (near) equality will emerge for V{D — > KK + tt's) vs. T{D — > 
tt's). Testing this expectation requires the measurement of final states with 
neutrals. 

Experimental studies of multi-body decays with more than one neutral meson 
in the final state (in particular, Dalitz plot analyses) also allow us to have dif- 
ferent handles on the studies of direct CP- violation in D-decays in and beyond 
the Standard Model [96] as well as on the dynamics of hadronic resonances 
governing these transitions (see, e.g. the E791 analysis [97,98]). 

6.3.2 D„ D+ ^ n+u, T+i/ 

The primary goal behind measuring leptonic decays, £+z/, or the Cabibbo 

suppressed versions, D'^ — > i'^u, with £ — /i, t, is the desire to extract the 
decay constants and fo. These quantities are important probes of heavy 
meson wave functions. In addition, these decay constants have been extracted 
from Monte Carlo simulations of QCD on the lattice with estimated uncer- 
tainties of about 20 percent on their absolute values and about 10 percent on 
their ratio. Improvements are expected for future lattice calculations. For the 
proper evaluation of these calculations, one wants to calibrate them against 
experimental results of similar accuracy. 

Currently, the branching ratios for Dg — > iv transitions have been measured 
by the CLEO collaboration with large uncertainties, Br{Ds — > ^iv) — 4.0l2;o ^ 
10~^ and Br{Ds tv) = (7±4) x 10^^. No measurement is currently available 
for other D mesons although there is an upper bound: Br{D^ fiu) < 
7.2 X 10~^. This can be explained by the A = 0.2 CKM suppression factor for 
the leptonic decays relative to those of the Dg. 
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Once the absolute vahics of fo or fn^ are known experimentally with about 
% accuracy or better then one will be able to feel more confident about ex- 
trapolating to the decay constants in the B system, Jb and Jb^, which are 
crucial quantities for a quantitative understanding of — B^ oscillations and 
the extraction of Vtd from them. 

Observing and measuring these transitions has always represented a highly 
nontrivial experimental challenge (and much more so for — > i'^i'), so the 
potentially exceptional performance for observing 1-prong D decays at i/MC's 
could allow them to make a significant contribution here even down the line. 
As a secondary goal one might even perform a detailed comparison of the 
rates for D ^ fiu and D ^ tu as a probe for new physics in the form of a 
non-minimal Higgs sector, for charged Higgs exchanges would affect the latter 
much more than the former. 

6.3.3 Inclusive Charm Hadron Decays 

Heavy quark expansions (HQE) allow the treatment of inclusive heavy flavor 
decays, including their non-perturbative aspects [99-103]. In addition to to- 
tal decay widths, other central quantities are inclusive semileptonic branching 
ratios and decay spectra for the different meson and baryon species. These 
techniques provide the basis for some of the most reliable methods for ex- 
tracting \Vcb\ and \\ub\ in B decays. Obviously one wants to cross check these 
methods in a system where the CKM parameters are known, namely the charm 
system, by testing how precisely |\4s| and \ Vcd\ can be extracted from semilep- 
tonic charm decays. In addition one can extract the size of the matrix elements 
of four-fermion operators that are of direct relevance in beauty decays and at 
the same time provide important calibration points for lattice simulations of 
QCD. 

No data of sufficient detail are available. The B factories (CLEO, BaBar and 
Belle) will significantly improve the situation, but might not achieve the de- 
sired experimental accuracy. Furthermore, it turns out that comparing neu- 
trino with charged lepton spectra in semileptonic decays provides us with 
particularly probing insights. 

One has to keep the following in mind. Since the expansion parameter is 
fJ'had/iTT'c with ^had ~ 0.7 — 1 GcV, ouc has to allow for uncalculated higher 
order contributions to modify the results significantly in charm decays. To 
have a handle on this complication, one needs to be able to perform detailed 
comparisons of the lepton spectra separately in Dg and Ac decays, 

which should be possible at uMCs. 

FOCUS and SELEX data will presumably yield precise lifetimes for 
baryons, but quite possibly not for the Qc- The latter is presumably the short- 
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est lived hadron in the single charm sector, with r(f2c) < 10^^^ sec; due to 
its different spin structure its lifetime is affected by different matrix elements 
than for the other baryons. It is also quite unclear whether next generation ex- 
periments like LHC-B and BTeV can measure such a short hfetime with good 
accuracy. A i/MC thus could make a relevant measurement that would serve 
as an a posteriori calibration of some theoretical tools. Furthermore, a whole 
new spectroscopy could be entered into, namely that of baryons carrying two 
units of charm: [ccq] . 

In principle, radiative inclusive (and exclusive) decays can also be studied. The 
predicted branching ratio for the short-distance contribution is tiny, Br{c 
wy) = (4.2 — 7.9) x 10"^^ [104], although two-loop QCD corrections could 
bring it up to 5 x 10^*^ [105]. This could have made it a sensitive probe of new 
physics as these processes occur in the Standard Model only at one loop. Unfor- 
tunately, the problem is that the long-distance effects can actually completely 
dominate this decay, enhancing it up to ~ 10~^, and these enhancements 
cannot be estimated model-independently. 

6.4 Searches for New Physics in Charm Decays 
6.4.1 D° - 5° Oscillations 

The phenomenon of meson-antimeson mixing has been studied both experi- 
mentally and theoretically for a long time as it provides an extremely sensitive 
test of the Standard Model as well as its various possible extensions. This is 
especially true for — mixing, as was already indicated. 

To study such oscillations one must tag separately the flavor of the produced 
meson and of the decaying meson. The charge of the primary lepton from a 
CC interaction uniquely tags the production sign of the charm quark method 
at an uMC This should easily be the cleanest and most efficient tag. It can 
be checked by a more conventional alternative method involving production 
of the charged D* mesons and studies of the decay chain D*^ D^{D^)7r^ 
[106,107], where anti-correlation studies of the charge of tt and decay products 
of D would reveal whether mixing took place. 

For charmed mesons with tagged production flavor, — oscillations 
are most cleanly probed through 'wrong-sign' semileptonic decays with the 
branching ratios: 



for Am/j = m(D°) — m^o, the D width and AFd the difference in the D° 




(111) 
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and mass widths. 

In principle, one can determine tiie flavor of tlie flnal state through charged 
kaons; mis-tags that happen due to doubly Cabibbo suppressed decays can be 
eliminated using a time-dependent analysis, as discussed below. 

The most recent experimental limits, which are from fixed target experiments 
at FNAL and from CLEO at CESR and combine tagging through 'wrong' sign 
leptons and kaons, read: 

rz) < 5x 10-\ 95%C.L.; CLEO[108] (112) 

-0.04 <yD< 0.06 , 90% C.L.; E791[109] (113) 

-0.058 <y'D< 0.01 , 95% C.L.; CLEO[108] (114) 

where 



Vd = VdCOsSkw -XDsmSxTT, (115) 

with Sktt denoting the strong phase shift between K'^tt" and 5° — > 

K'^7r~ (see [110] for the recent analysis), and 

ycp = 0.0342 ±0.0139 ±0.0074; F0CUS[111] (116) 

where ycp — Vd in the Standard Model. Since possible new physics effects 
or hadronic uncertainties will affect these experiments differently, a careful 
analysis to extract the true values of Arrio and AF/j from the data should be 
performed [112]. The B factories at Cornell, SLAC and KEK will refine the 
search for — D° oscillations to an expected sensitivity of ro ~ fewxlO"^ 
[113]. 

While the Standard Model undoubtedly predicts slow — oscillations - 
xd, Vd ^ i - there is considerable uncertainty in the numerical predictions. 
A conservative Standard Model bound is given by [114,115] 

rolsM < 10-" ~ Vd, xn\sM < 10"" . (117) 

Bolder predictions have been made that xd and yo cannot exceed 10"^ [114,115] 
and therefore < 10~^ within the Standard Model. On the other hand, new 
physics could enhance xd up to, and actually even above, the present bound, 

xd\np ~ 0.1 , (118) 

without violating any other limit and while leaving yo unaffected. Examples 
of such new physics processes include various supersymmetric models [116] 
(including SUSY models with quark-squark alignment that actually require 
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Ami) close to the current experimental bound) [117], models with singlet up 
quarks [118], various leptoquark models [119], and multiscalar models with 
[120], and without [121], natural flavor conservation. Any experimental effort 
to lower the current limit on Am^) is essential in determining the available 
parameter space for many possible extensions of the Standard Model! 

The cleanest way to probe for D'^ — oscillations is to analyze the time 
evolution of transitions into 'wrong-sign' leptons: 

T{D%t) ^ tX) oc e-'/^'^xl (-^y . (119) 

Here we have invoked the AQ — —AC rule of the Standard Model which 
makes oscillations the only source for wrong-sign leptons. 

Since one is embarking on a search for new physics, one should generalize 
equation 119 to allow for a violation of the AQ — —AC rule, giving: 

r(D°(i) ^£+X) oce-*/^^x 



(1 \ X 1 ^ p 

1 + -ATntj \Pwrong\^ + - {AmDtY - -AVDtRe-pwTong + Am nthn- Pwrong 



(120) 



where 



denotes the ratio of AC = AQ to AC = —AQ amplitudes, 

\Di,2)=p\D')±q\D') (122) 

relates mass and flavor eigenstates, and the oscillating functions multiplying 
the usual e^*/^° term have been expanded in powers of the proper time t since 
^D, Ud ^ 1- The AC = AQ term has no t dependence beyond that of e~*/'^^, 
the pure oscillation term has a dependence, while the interference between 
the two generates a term linear in t. 



The violation of the AQ = —AC rule arises even within the Standard Model 
for the decays D° — > X+tt" due to doubly Cabibbo suppressed transitions 
(DCST) producing the direct decay — > X+tt", with a branching ratio 
Br{D^ K+TT-) = (2.8 ± 0.9) ± 10"^ [110,122], and thus mimicking the 
signal for DD mixing. The equation corresponding to Eq. 120 is: 

r{D\t) K+n-) oc e-^o°*tan%c|pi^^r 
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X 



i + Iaf t (^^^^)' ^^^^ 

2 4tan4ec|p;f.P 2tan2ec|p;f.| 



Amot 



where 



tan^6'c|pK^| 



T(L'0 ^ X+TT-) 



(123) 



(124) 



is the fraction of wrong-sign decays. 



One can also search for hfetime differences in certain well-chosen decay 
channels in order to probe the contributions to oscillations from the term 
of Eq. 111. With CP invariance holding (at least) to good approximation, CP 
eigenstates can be treated as mass eigenstates. While D*^ — > K~^K~, tt+tt" 
will then exhibit r+, D° ^ Ks4>, Kslo, Ksp, KsT] etc. will be controlled by 
r_, where r+ [r_] denotes the width for the CP even [odd] state and 



Ar 



(125) 



Furthermore the width for D° K tt"^ is approximately given by (r+-|-r_) /2 
[111,112]. 



6.4-2 CP Violation in D Decays 

There is a wide field of potential CP violation in D decays that can be dis- 
cussed in close quahtative analogy to B decays. 

CP asymmetries that necessarily involve — D^ oscillations can arise in final 
states that are CP eigenstates, like K^K~ or 7r"'"7r~: 

r(D°(t) ^ K+K-) oc e-^°* ( 1 + sinAmDi • Im^p/^+x- j 



P . / Am,i)t t q \ , ^ 

~ e-^-* 1 + — ^ Im^PK+K- . 126 

\ Td td p J 



With xd\sm < 10 ^ and Iuij^pk+k-Ikm ~ 0(10 '^), one arrives at an asym- 
metry of only around 10~^, which would likely be too small to measure even 
at a z/MC. Yet with new physics one conceivably has xd\np < 0.1 and 
lm^PK+K-\NP ~ ^'(lO^^), leading to an asymmetry that could be as large as 
of order 1%. 
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Likewise, one can search for CP violation by comparing the proper time distri- 
bution of Eq. 123 for the doubly Cabibbo suppressed transitions K^t:~ 
with that for decays: 

T{D^{t) K-7r+) oc e-^o°*tan^^c|Pi^^|^ 



X 



2 ^ Atan^ec\pKn\'^ 2tan2 KpIPkJ 



tan^ 9c\pKn\ \P\PKn\ 

where 



(127) 



ta..-.c-fe, 4g:^:::j . (i28) 

In such new physics scenarios one would expect a considerably enhanced asym- 
metry - perhaps as large as 1%/ tan^ 9c ~ 20% - but at the cost of smaller 
statistics. Hoping for an asymmetry of several percent is more realistic, though. 
Effects of that size would unequivocally signal the intervention of new physics! 
One should note that these rough estimates are based on xd — 10~^ which 
would correspond torn — 10"^. This implies that, even if oscillations have not 
been found on the = 10~^ level in semileptonic decays, a CP asymme- 
try of several percent (or conceivably ten percent) could still be encountered 
in D° K+7r-\ 

Direct CP violation can occur as well. There are actually two types of effects: 
differences between partial rates for CP conjugate transitions 

_ T(D ^ /) - T(D ^ /) 

^'■^ - r(c^/) + r(B^/) ■ 'i'") 

and asymmetries in final state distributions such as, e.g., Dalitz plot popula- 
tions. 

Strong final state interactions play an important part in both cases: in the 
former they must induce the phase shifts that are essential to make a differ- 
ence observable; in the latter they can very significantly affect the observable 
asymmetry. The existence of resonances in the neighborhood of the charmed 
meson mass is proof that hadron dynamics is active in this energy region and 
will affect the weak decays of charmed particles. 

The good news is that whenever there are CP violating weak phases one 
can count on final state interactions to make them observable. The bad news 
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is that interpreting a signal as evidence for new physics will pose a highly 
nontrivial theoretical challenge. 

The Standard Model with the CKM ansatz can induce direct CP asymmetries 
only in Cabibbo suppressed channels. Model-dependent estimates usually pre- 
dict direct CP asymmetries to be of the order of 10~^ but, exceptionally, they 
could reach the 10~^ level [123]. A measurement of the branching ratios for 
all related channels - in particular also those with neutral hadrons in the final 
state, as sketched above - would enable us to constrain the strong phase shifts 
quite significantly. A i/MC will have a significant advantage in this respect! 



6.4-3 T Odd Correlations in Ac Decays 

One special feature of uMCs is represented by the production of Ac in CC and 
NC reactions: 

uN uAcX or iiAcX. (130) 

This allows novel studies of various Ac decay form factors with a Q"^ range 
extending to well above mj^^. Yet even more intriguing and promising would 
be a detailed analysis of the final state in its semileptonic decays: 

A+ ^ i+ueA. (131) 

With the parent c quark being left-handed one expects the Ac to emerge in a 
highly polarized state. The usual valence quark description actually suggests 
that the Ac polarization is completely carried by its c quark; i.e., a left-handed 
c quark fragments into a left-handed Ac. Yet even with unpolarized Ac one can 
form an experimentally observable T odd correlation 

Ct+- = (<JA • (pa X Pi)) (132) 

connecting spin and momentum of the daughter hyperon with the lepton mo- 
mentum. 

In a general experimental process, observing a non-vanishing value for a T 
odd correlation does not automatically establish that T (and CP) invariance 
is violated since in general final state interactions could fake such an effect. 
However, this problem does not occur for Eq. 131 since it cannot be affected 
by either strong or electromagnetic final state interactions! This is analogous 
to the well known situation in ji^vii^ ( vs. Kl fi'^uii^). Like there, 

the CKM ansatz cannot generate an observable effect here, yet certain new 
physics scenarios can. An effect of order 10"^ is not inconceivable, particularly 
if the channel Ac — > ruA could be studied [124]. 
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6.5 Summary on Charm Decay Physics at vMCs 

The research program at a uMC is hkely to improve our knowledge and un- 
derstanding of charm decays quite significantly even ten years from now: 

• It would fill out many white spots on our map of Ac, 5c and VLc decays 
by measuring many new relative and absolute branching ratios, including 
for final states with more than one neutral particle. 

• It would allow the measurement of , /lu, tu in a very clean 
environment. 

• It could probe for — oscillations and for CP asymmetries involving 
them with superbly clean systematics. It would significantly improve on the 
sensitivity that can be obtained at B factories for such phenomena. 

• It would enable us to search for direct CP asymmetries in many different 
channels and at the same time provide us with information that could help 
us in properly interpreting a signal. 



7 Neutrino Oscillation Experiments with a Muon Storage Ring/Neutrino 
Factory 

7. 1 Status of Neutrino Oscillations at the Time of uMCs 

In a modern theoretical context, one generally expects nonzero neutrino masses 
and associated lepton mixing [9]. Experimentally, there has been accumulating 
evidence for such masses and mixing. All solar neutrino experiments (Home- 
stake, Kamiokande, SuperKamiokande, SAGE, and GALLEX) show a signifi- 
cant deficit in the neutrino fiuxes coming from the Sun [125]. This deficit can be 
explained by oscillations of the z/g's into one or more other weak eigenstates, 
with Am^yj of the order 10"^ eV^/c^ for solutions involving the Mikheev- 
Smirnov-Wolfenstein (MSW) resonant matter oscillations [126,127] or of the 
order of 10~^° eV^/c^ for vacuum oscillations. Accounting for the data with 
vacuum oscillations (VO) requires almost maximal mixing. The MSW solu- 
tions include one for small mixing angles (SMA) and one with essentially 
maximal mixing (LMA). 

Another piece of evidence for neutrino oscillations is the atmospheric neutrino 
anomaly, observed by Kamiokande [128], 1MB [129], SuperKamiokande [130] 
with the highest statistics, and by Soudan [131] and MACRO [132]. This data 
can be fit by the inference of oscillations with Am^^^ ~ 3.5 x 10^'^ 

eV^/c^ [130] s-iid maximal mixing siv?26atm = 1- The identification = Uj- 
is preferred over i/^ = ^sterile about the 2.5(7 level [133], and the identifica- 
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tion fx = Ve is excluded by both the SuperKamiokande data and the Chooz 
experiment [134,135]. 

In addition, the LSND experiment [136] has reported observing — > Pg and 
i^/i — ^ oscillations with /^^^mj^gj^j^ ~ 0.1 — 1 eV^/c'^ and a range of possible 
mixing angles, depending on Am^gj^jj. This result is not confirmed, but also 
not completely ruled out, by a similar experiment, KARMEN [137]. Inclusion 
of the signal reported by LSND with the other two pieces of evidence would 
imply three distinct mass differences and hence four neutrinos. Some proposals 
for the form of the mixing matrix invoke only 3 generations of neutrinos to 
account for all signatures, while others invoke a fourth sterile neutrino. 

A number of fits have been made to the existing neutrino data. The fit by 
the SuperKamiokande collaboration to its data yields a minimum in the at 
sin^(2^at^) = 1, with an allowed region of 0.8 < sm'^{29atm) 1- In terms of 
the basic angles in the Icpton mixing matrix, this implies that 623 is close to 
71 /A and allows a small, nonzero ^^13, consistent with the bound from CHOOZ. 
As will be discussed below, a major physics capability of the muon storage 
ring/neutrino factory is the ability to measure ^13. 

There are currently intense efforts to confirm and extend the evidence for neu- 
trino oscillations in all of the various sectors - solar, atmospheric and accelera- 
tor. Some of these experiments are running; in addition to SuperKamiokande 
and Soudan-2, these include the Sudbury Neutrino Observatory, SNO, and 
the K2K long baseline experiment between KEK and Kamioka. Others are in 
development and testing phases, such as BOONE, MINOS, the CERN-Gran 
Sasso program, KAMLAND, and Borexino [138]. Among the long baseline 
neutrino oscillation experiments, the approximate distances are L ~ 250 km 
for K2K, 730 km for both MINOS, from Fermilab to Soudan and the proposed 
CERN-Gran Sasso experiments. The sensitivity of these experiments is pro- 
jected to reach down roughly to the level Am^ ~ 10~^eV^/c"'. Experiments 
that are planned as part of this program include [139] ICANOE and OPERA 
[140]. Although they are expected to begin operation after MINOS, they will 
involve somewhat different detector designs and plan to focus on establishing 
T appearance. This, then, is the program of research for the next several years. 



7.2 Oscillation Experiments at vMCs 

Although a neutrino factory based on a muon storage ring will turn on several 
years after this near-term period in which K2K, MINOS, and the CERN-Gran 
Sasso experiments will run, we believe that it has a valuable role to play, given 
the very high- intensity neutrino beams of fixed fiavor-pure content, including, 
in particular, and Pg beams as well as the conventional i/^ and beams. 
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The potential of the neutrino beams from a muon storage ring is that, in 
contrast to a conventional neutrino beam tt^ / decay, is primarily with 
some admixture of z/g's and other flavors from K decays, the neutrino beams 
from the muon storage ring would high extremely high purity: /x^ beams would 
yield 50 % and 50 % Pe, and so forth for the charge conjugate case of ji^ 
beams. Furthermore, these could be produced with extremely high intensities, 
of order 10^° to 10^^ neutrinos per year. 

Given the form of the oscillation probabilities, a neutrino beam for an oscilla- 
tion experiment would optimally be made from the lower end of the spectrum 
of energies considered in this report. Because of the lower requirements on 
beam focusing and acceleration, making a muon storage ring for the purposes 
of an oscillation experiment has been proposed as a first step towards develop- 
ing a muon collider. Energies being considered range from 20 GeV to 50 GeV, 
and the geometry of the final muon ring is very different from a traditional 
collider ring, in that the straight section that points to a neutrino experiment 
comprises between 25 and 40% of the "circumference" of the ring. 

The types of neutrino oscillations that can be searched for with the neutrino 
factory based on a muon storage ring, along with the final state charged lepton 
species for neutrino-nucleon DIS that tags the interacting neutrino fiavor, are 
listed below for the case of the i/^^e beam from /x~, decaying as ii~ —>■ i/^e~Pe: 

> li~ (survival); 

e~ (appearance); 

> T~ , T~ (e~,/x~)... (appearance*); 
e (survival); 

> (appearance); 

■ T+; (e+,/x+)... (appearance*); 

where the asterisks denote that the tau appearance signatures may be some- 
what indirect in involving detection of the leptonic daughters of tau decays 
rather than the decay vertices of the taus themselves. 

It is clear from the list of processes above that, since the beam contains both 
neutrinos and antineutrinos, the only way to determine what the parent neu- 
trino was is to measure the charge of the final state lepton. The — > 
oscillation will produce a wrong-sign e~ as will the z/^ Vr oscillation fol- 
lowed by r decay to e~. The easiest wrong-sign lepton signatures to detect 
arise from the oscillations u^. — > P^, giving a /x"*", and from u^. — > u^-, giving 
a r"*" which will decay part of the time to If one is searching for r final 
states, muon storage ring energies above 30 GeV should be used to minimize 
the threshold kinematic suppression. 

To get a rough idea of how the sensitivity of an oscillation experiment would 
scale with energy and baseline length, recall that the event rate in the absence 
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of oscillations is simply the neutrino flux times the cross section. First of 
all, neutrino cross sections in the region above about 10 GeV (and slightly 
higher for r production) grow linearly with the neutrino energy. Secondly, the 
beam divergence is a function of the initial muon storage ring energy; this 
divergence yields a flux, as a function of Od, the angle of deviation from the 
forward direction, that goes like 1/^^ ~ E'^. Combining this with the linear E 
dependence of the neutrino cross section and the overall 1/L^ dependence of 
the flux far from the production region, one flnds that the event rate goes like 

To set the scale, consider an experiment that sees 2x10^° 30 GeV muon decays 
in a straight section pointed at a detector 2800km away. In the absence of 
oscillations, the Vf^iVe) charged current rate would be 52,500 (22,600) events 
per 10 kton [7]. Figure 16 shows the relative Ve and statistics for two 
conflgurations: one is a 20 GeV /x"*" ring with a detector at 2800 km, the other 
is a 50 GeV ring with a detector at 9100 km, which is close to the distance 
from either Fermilab or CERN to Kamiokande. 

Now recall the general formula for the probability of a two-species neutrino 
oscillation in vacuum, say of —>■ i/^: 

P{u, u,) = W\U,s\'sm' (^%^) (134) 
= sm'{29^s) sm\9,s) sin' l^^^) ■ 



Where the sin^ term is small, it can be expanded to give a factor of {L/EY 
that cancels the L dependence in equation 133 and reduces the E dependence 
to a linear factor. Underlying these considerations of optimal energy is the fact 
that, even if one designs an experiment for r appearance, the overall event rate 
of detected r's may be rather small. 

It is quite likely that, by the time a neutrino factory turns on, Am^^^ and 
sin^(2^23) will be known at the 10-30% level. Although a neutrino factory could 
undoubtedly improve the precision on those two parameters, the novel physics 
that can be addressed is a determination of 6*13, and the sign of Am^^^ . By 
using matter effects, and a comparison of Va versus Va oscillations (by switching 
the muon storage ring from //~ to //"*"), the sign of Am^^^ can be determined. 
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Fig. 16. The upper plot shows fe and z7^ event rates per GeV as a function of E/L 
for two different experiments. The lower plot shows the L/E dependence of the 
oscillation probability, assuming the largest 5m? is 2, 3. 5, or 5 xlO^'^ eV^/c^. 

7.3 Matter Effects 



With the advent of the muon storage ring, the distances one can place de- 
tectors will become large enough that, for the first time, matter effects can 
be exploited in accelerator-based oscillation experiments. Simply put, matter 
effects are the matter-induced oscillations which neutrinos undergo along their 
flight path through the Earth from the source to the detector. Given the typ- 
ical density of the earth, matter effects are important for the neutrino energy 
range E ~ 0(10) GeV and Am^^^ ~ 10^'^ eV^/c^ values that are relevant for 
z/MC long baseline experiments. 

Follow-up studies to initial discussion [126] of matter-induced resonant neu- 
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trino oscillations include an early study of these effects that assTimed three 
neutrino generations [141] and a discussion of the sensitivity of an atmospheric 
neutrino experiment to small Am? due to the long baselines and the necessity 
of taking into account matter effects was discussed[142]. Many analyses were 
performed in the 1980's of the effects of resonant neutrino oscillations [127] 
on the solar neutrino flux, and matter effects in the Earth [143,144], and of 
matter effect on atmospheric neutrinos [145]. Work continues[146,147] on mat- 
ter effects relevant to atmospheric neutrinos. Early studies of matter effects 
on long baseline neutrino oscillation experiments [148] have been extended to 
cover neutrino factories [7, 149], [150,151] 

In recent papers by one of the authors (RS) and 1. Mocioiu, calculations were 
presented of the matter effect for parameters relevant to the possible long 
baseline neutrino experiments envisioned for the muon storage ring/neutrino 
factory [152,153]. In particular, these authors compared the results obtained 
with constant matter density along the neutrino path versus results obtained 
by incorporating the actual density profiles. They studied the dependence 
of the oscillation signal on both E/ Aml^^ and on the angles in the leptonic 
mixing matrix, and commented on the influence of Am^o; and CP violation on 
the oscillations. Additional recent studies are listed in the references [154,155]. 



7.4 Detector Considerations 

In order to measure oscillation parameters that describe the transitions above, 
one would ideally want a detector that could identify the existence and flavor 
of any outgoing lepton from the neutrino interaction, as well as the hadronic 
and leptonic energy in the event. When measuring very small oscillation prob- 
abilities, however, backgrounds must be taken into account. Naively one would 
think that simply detecting a muon of opposite charge to that in the storage 
ring is a signal for the electron neutrino oscillating. However, pions and kaons 
are produced copiously in neutral and charged current neutrino interactions, 
and if one decays to a muon before it interacts in the detector this can con- 
stitute a significant background. At higher energies charmed mesons are also 
produced which decay immediately to muons 10% of the time. Ultimately, 
detectors will need to be designed that have sufficient resolution on both the 
energies and the angles of the final state lepton and hadronic shower to be 
able to remove these backgrounds. 

Although detectors exist that could identify all of the final state leptons and 
their charges, the challenge is to make them on the several kiloton scale. If 
the largest Am^ is in the LSND region, there will undoubtedly be more work 
done to optimize relatively low mass detectors that emphasize tau appearance. 
However, the only detectors which have thus far been proposed on the 10-40 
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kiloton scale are for detecting wrong sign muon events. The two detector tech- 
nologies that have been considered in detail for oscillation experiments for a 
muon storage ring will now be discussed in turn: one is a magnetized sampling 
calorimeter such as the one used by MINOS [156] and the other is a liquid 
argon time projection chamber (TPC) combined with a muon spectrometer 
such as the one proposed by the ICANOE [157] collaboration for the CERN 
to Gran Sasso neutrino beam. 

7.4-1 Magnetized Sampling Calorimeters 

Magnetized sampling calorimeters consist simply of alternating layers of mag- 
netized steel and readout, where the readout traditionally consists of scintil- 
lator and/or drift chambers. 

The charged particle efficiency of the readout planes can be close to 100% so 
the performance of the sampling calorimeter depends primarily on the sam- 
pling frequency of the detector. A steel/scintillator sandwich with sampling 
every 5 cm of steel would have a fractional hadron energy [H) resolution of 
approximately [156] 



With fine enough transverse segmentation, the hadron angular resolution is 
dominated by the hadron energy resolution. The muon energy and angular 
resolution are expected to be much better than for the hadronic shower. 

Although separating z/g charged current events from neutral current events 
is difficult and determining the charge of the outgoing electron impossible in 
this detector, a muon in the final state can be easily and efficiently detected, 
and its charge, momentum, and initial outgoing angle can be determined once 
the muon traverses enough steel to be spatially separated from the hadronic 
shower. Kinematic cuts can be made on the muon momentum and its compo- 
nent transverse to the hadronic shower to reduce the background from charm 
production. With signal efficiencies from 25 to 30%, the backgrounds can be 
reduced to a level of 10~^ to 10~^, depending on the neutrino energy. At higher 
energies the backgrounds are larger but the faster improvement in the back- 
ground rejection actually causes a reduction in the background contributions 
to the analysis [158]. 

74.2 Liquid Argon TPCs 

The ICANOE-type detector would consist of a large volume of liquid argon in- 
strumented with time projection chambers (TPCs), followed by a much thin- 
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ner volume of magnetized steel where a muon's charge and momentum can be 
determined. The TPC would have very small wire spacing (3mm) and would 
act much like an electronic bubble chamber. Electron neutrino charged current 
interactions could be distinguished from neutral current interactions, although 
the electron charge could not be measured. By breaking up the event samples 
into four distinct classes - right sign muons, wrong sign muons, electron-like 
events and neutral current events - one could fit all four distributions simul- 
taneously to determine oscillation parameters, u^s might also be identified 
on a statistical basis by looking at the acoplanarity distribution in the event 
sample. 

The energy and angular resolutions of all the final state particles would be 
extremely good, e.g.. 



for the hadronic energy and 150 mrad for hadron shower angles. However, the 
ability of this detector to see wrong-sign muons would depend primarily on 

the segmentation between the liquid argon and the magnetized spectrome- 
ter and on the thickness of the spectrometer itself. The thinner one makes the 
spectrometer, the more likely one is to have backgrounds from charge misiden- 
tification. The thicker the spectrometer, the less room there is for the liquid 
argon in a given volume. The thinner the liquid argon, the higher the accep- 
tance for low energy muons (since muons lose approximately 210 MeV/m in 
liquid argon [159]), but the less target volume one has overall. Clearly opti- 
mization of this geometry is needed, and will depend somewhat on the energy 
of the muon storage ring. 

7.4-3 Muon Detector Conclusions 

Although the two types of detectors have different strengths, detailed stud- 
ies [160,149,150,153,161,3,4] have shown that both either would be adequate 
to make precise measurements of |(5m^|, sign (Sm^), and sin^ 623; and to extend 
the sensitivity of sin^ 6*13 by 1-2 orders of magnitude in the scenario where the 
largest Am^ is described by the atmospheric neutrino anomaly. 

7.4-4 Tau and Electron Detectors 

Alternate technologies must be employed to achieve electron or tau identi- 
fication event-by- event, or electron or tau charge measurements. If LSND is 
confirmed and the largest Am^ would suggest baselines on the order of tens 
of km, then a much higher premium will be placed on designing detectors 
that can do tau and electron charge determination, and they will not have 
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to be as massive. At these short baselines, detectors on the 1 kiloton scale 
could be quite adequate to make precision measurements on z/g Uj. and 
i^ii ~^ ^T- Even if LSND is not confirmed, efforts to make massive tau and 
electron charge identification detectors should not be abandoned since these 
two channels still comprise a large part of the mixing matrix and should be 
researched to confirm our understanding of neutrino mixing. 

One category of new detectors uses thin (~ 100 //m) sheets of emulsion com- 
bined with low-density (~ 300 \im) spacers, and thin sheets of metal to give 
the detector mass. With emulsion one can measure the kink that occurs when 
a tau decays by comparing the slope of a track before and after the spacer. 
Such a geometry, with lead as the mass, is described in reference [162]. This 
would be very useful for identifying taus and electrons. However, for charge 
identification one needs to introduce a magnetic field. This could be done 
using an extremely large external magnet, such as the one used in ATLAS, 
and thin steel plates, or by using a coil and magnetized steel to make the 
mass [163]. Since the overhead for analysis of each event is high in this sort of 
detector, one would place it in a region where the tau appearance probability 
is maximized. 



7.5 Conclusions on Neutrino Oscillation Studies at vMCs 



In conclusion, neutrino masses and mixing are generic theoretical expecta- 
tions. The seesaw mechanism naturally yields light neutrinos, although its 
detailed predictions are model-dependent and may require a lower mass scale 
than the GUT mass scale. Current atmospheric neutrino data is consistent 
with maximal mixing in the relevant channel, which at present is favored to 
be — ^ Vr- Even after the near-term program of experiments by K2K, MI- 
NOS, the CERN-Gran Sasso experiments, and mini-BOONE, a high-intensity 
neutrino factory generating 10^° — 10^^ neutrinos per year will add greatly to 
our knowledge of the neutrino masses and mixing matrix. Ideally, the muon 
storage ring should be coupled with two long-baseline neutrino oscillation ex- 
periments, located at different baselines, that can take advantage of matter 
effects to amplify certain transitions and with a massive detector that will 
identify /x's and r's with charge discrimination. In particular, it should be 
able to measure Am^^^ and sin^(2^23) to the level of several per cent and also 
give important information about the sign of Am^^^ and about sin^(2^i3). 
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8 Summciry 



Beams from uMCs have the potential to provide vast improvements over to- 
day's conventional neutrino beams from tt/K decays. They are much more 
intense, have a much smaller transverse extent and produce precisely pre- 
dictable beam spectra. 

Oscillation experiments may extend even to intercontinental baselines, while 
unprecedented event statistics approaching of order 10^ to 10^'^ precisely con- 
structed DIS events will open new regimes of neutrino interaction physics. The 
extraordinary rates will enable use of active vertexing targets surrounding by 
a high resolution spectrometer and calorimeter. The first polarized targets for 
neutrino scattering can be substituted for special studies. 

Highlights of uMC physics program include: 

• substantially extending the reach of accelerator-based experiments to study 
neutrino oscillations; 

• measurements of the CKM quark mixing matrix elements \Vcd\, \Vub\, \Vcs\, 
and \ Vcb\ in inclusive high scattering, with few percent accuracies achiev- 
able for the first two; 

• a realistic opportunity to determine the detailed quark-by-quark structure 
of the nucleon; 

• mapping out the quark-by- quark spin structure with polarized targets, and, 
perhaps, determining the gluon contribution to the nucleon's spin; 

• some of the most precise measurements and tests of perturbative QCD; 

• tests of the electroweak theory through measurements of sin^ 9w with frac- 
tional uncertainties approaching 10~^; 

• a new realm to search for exotic physics processes; 

• a charm factory with unique and novel capabilities; 

• a new laboratory for the study of nuclear physics with neutrino beams. 

The potential experimental capabilities of i/MCs reach so far beyond present 
neutrino programs that physics surprises not touched upon in this report can 
be expected. To see, we must build the machines. 
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